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FLUCTUATION ACTION IN THE SIMPLEST PARAMETRIC SYSTEMS 


I. Tikhonov 


(Moscow) 


The statistical characteristics of the output, random signal of a parametric system 
are obtained for the case when the external action and the variable parameter are corre- 
lated, stationary, normal, random functions, The system under consideration is described | 
by a first-order, linear, differential equation, 


In many cases it is not possible to rigorously solve the problem of the action of random signals in non- 
linear systems, necessitating, as a consequence, the introduction of certain simplifications, At the same time 
it is sometimes possible to reduce the problem to the solution of a first-order, linear, differential equation, 
with a variable coefficient of the form 


+ 


where € (t) and n(t) are random functions, Henceforth, we shall assume that these functions are correlated, 
stationary and normal, with known characteristics. 


Let us introduce the following notations: 


m, = CE(t)>, me = (t) = <[E — my] [g(t + — 
= — m,] [E(t + t) — m,]>» 
(2) = — ma] [4 (¢ + *) —my]), 
where the skew brackets denote the mathematical expectation operation, i.e., the operation of averaging over 
many different realizations, 


As is well known [1], the general solution of Equation (1) can be expressed in the form of the sum of the 
general solution y, of a homogeneous equation and the particular solution y, of an inhomogeneous equation, 
For the initial conditions y = yp, t= ty, we can write 


Let us note that this solution of Equation (1) can be obtained in a different way, The solution of the 
differential equation 
d: 
ar + = exp[—an(t)}, 


where a is some parameter, for the initial conditions z = yg, t=t, has the form: 
t t t 


= + = das) + oxp|—an(u)— da] de (a 
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From a comparison of Expressions (2) and (3), it follows that 


This result permits a significant simplification in the calculation of statistical characteristics y(t), In 
particular, for the calculation of various moments it is sufficient to utilize the expression for the n-dimensional 
characteristic function of a normal random process: 


uml 


Here m, =m(t,), o, = o(t,) are, respectively, the mean and mean-square values of the random 
process at some instant of time t,,Ry, (r)=R(ty— ty) is the coefficient of correlation between two values 
of the random process at times t, and t,, 


Let us calculate the mean value and the correlation function of y(t), From Relationship (4) we have 


m (t) = <y (t)> = yo<exp(— >] 


Insofar as the linear transformations of a normal process yield a normal process too, then the first term 
of Equation (6) represents the value of the one-dimensional characteristic function of the normal process 


t tt 
= my (t— te) = (t, t) oth 


at the point i,e., 


§,(/) = Cexp dz,)> 


= exp [—m, (t— + — a) (7 


be 
Analogously, making use of Expression (5), we obtain 


t 
94 (ja, j) = [—an(u) (2s) = exp {— fam, + m, (¢—u)]+ 


+ + (t, u) (Ru (u — a) day + (t, u)}} 


Consequently, 


= 9, (ja, du = 


= \[—m+ (t, u) (u — 2) x 


x exp[—m,(t—u) +5 


1 
u 


R, 2) dx, . (8) 
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Substituting (7) and (8) into Equation (6), we obtain 
tt 


m (t) = yo exp|— my (t — + — 2) dy da] + 


+ — Ras (u — 24) dag] x 


x exp[—m(t—n) + Lot | (9) 


The correlation function of y(t) is equal to © 


K (ty, tg) —= May (try ta) — mm (ty) (ty) = <y (ta)> — m (ty) m (ts), (10) 
where Mg (t;, t,) is a two-dimensional moment, For it we obtain from Formula (4). the expression 


Ma (tas ta) = (ta) 21 —[ (ta) 22 (4)? |, — 


In the calculation of the separate summands of the right side, one should make use of the expressions for 
the three- dimensional and four-dimensional characteristic functions of a normal process, which can be derived 
from Equation (5) for n=3 and n=4, 


Making use of the equations obtained in this manner, it can be shown that 


t 
(ta) = yg exp {— + + di + 


+ ox (ty, to) o¢ (tes to) — 21) dos + dnd 
ti, t,t, 


t 
20 (ta) {— cm, — my (ty + ts — — + 


(12) 


+5 [ad + — anand + | — dy + 


e*e 


ts 
++ to) | Raa — 21) dary + 2) | — + 
t, u 


+ (ths tp) o¢ (ta, Ry — dy du. 
From here 


(ts) 20 Yo \[— ma + (tr (u—2,)dx, + 
t, ° 


ty 
+ 020% (tay u) | — 2%) day] exp {— my + ta — to + (ty fe) + 


(tay + oy (t1, bo) (te, Ry (% —%) da, dz,} du. (13) 


Interchanging t and t, at this point we obtain the third term of the right side of Equation (11), The last 
(fourth) term in this expression can be written in the following form which is more convenient for calculation: 


Applying Equation (5) with n = 4 and carrying outthe necessary computations, we obtain 


(t1) 22 (t2)> | = (u—v) + 


t 
+ | mz — 4) (U — dx, — (ty, Ra (u— dz | x 


x [m, — 1) (ta, — %) dz,|} x 


x exp {— + + $08 (ty a) + $8 (ter + 
t, 
+ (ty, (ty, Ra (ea der dr (a4) 


Substituting the obtained Expressions (12) - (14) into Equation (11) and then into (10), we obtain a general 
formula for the correlation function, Since this equation turns out to be very cumbersome, we will limit our- 


selves here to examining certain particular cases, 


1, Assume zero initial conditions (y) =0, tp=0) and €(t)=my, Then Equation (1) becomes a first-order, 
linear, differential equation with a random external force, In this case we have from Equation (9) 


t 


The mean value of m in a stationary condition (t—> oo) will be finite and equal to m,/m, only when 
m> 0, i,e., if the system ossesses positive damping. We can, therefore, say that the system being considered 
is stable relative to its mean value only when m, > 0. 


For a two-dimensional moment we obtain from Equation (14) 


+ m5] exp [— m, (t; + tg — u —v)] dudv. 


0 
Making use of Equations (10) and (15), it is easy to see that the correlation function is equal to 


K (ty, te) met exp my (ty + tg —u — v)] du dv. 


00 
Introducing new variables x u, y=tg— Vv, T =ty— tg, t= ty, we can write 


—t 


tt 
k(t,t—*)= \ Ry (t + y exp [— m, (z+ dxdy. 


Let Rg(r+y— x) =exp (— 6|1r+y—x| ), for example. Carrying out separate integrations for r > 0 
(Fig. 1, a) and r < 0 (Fig. 1, b), by splitting up the regions of integration, in each case. into two sub-regions, 
depending on the sign of the argument of the coefficient of correlation, we obtain 


u 
uv 
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t+ bem —4)}, =>0. 
From this equation it follows that when my > 0 there is a stationary value of the correlation function 


2, Let us suppose n (t) = 0. In this case the mean value will be, according to Equation (9) 


Substituting in the variables xg— xy=s, xX, =x, and carrying out the integration, we obtain 


t—t, 
m (t) = yy exp[—m, (t— to) +o} (t — ty — 8) Ry(s) ds}. 


To determine the stationary condition we must let (t— t)—*» oo, At the same time, it is apparent from 
the equation that the mean value for the stationary con- 


dition is equal to zero only if 
m, — | Ry (s)dt > 0. ae) 
0 
For the reverse inequality the mean value tends to 


infinity and the system becomes unstable in this sense. 


Making use of Equation (12) it is easy to see that, 
for the given conditions, the correlation function is equal to 


(tay ta) = (ta) (ts) foxp[2og te) (tay R, (2% —%) dz, dzq|—1} - 


tt, 
From this we obtain the dispersion 


o* (t) = ygm* (t) {exp [223 (% — dx, dz) | —1} = 
= yim? (t) {exp [ 802 (t—t,—s)R, (8) ds) ‘ 


Insofar as for finite values of my and oy, for sufficiently large t > 0, the inequality 
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[ims — of | (s) ds} (¢— < (8) ds) 
0 0 
is always satisfied then, beginning at some instant of time, the dispersion will grow without limit and the system 
being examined (for zero initial conditions) turns out to be unstable with respect to the dispersion, despite the 
fact that when Condition (16) is satisfied the mean value tends to zero, 


There is a physical explanation for such a somewhat unusual result, Since the fluctuations € (t) are assumed 
to be normal, with mean value m and dispersion of, then there exist random intervals of time Aj (Fig. 2) 
during which €(t) > 0, and y(t), subject to the initial deviation 
ay Yo, dies away, However, besides the intervals A‘; there are inter- 
vals A‘, during which €(t) < 0, and, therefore, during them the 


a deviation y(t) grows, Despite the fact that the general time £ A}, 


ro ay during which €(t)> 0, is greater than the general time £A";, 
“Say ¢ during which €(t) < 0, the dispersion of y (t) grows with time as 
r a result of the fact that for €(t) w 0 the rate of growth of y(t) is 
ig. 2 
very large. 


3, Let the initial conditions be zero (yy = 0, ty = 0) and the random functions £ (t) and  (t) be noncorre- 
lated, Then from Equation (9), we have . 


m() =m, \oxp [— (¢—u) + d= 


= My my + (v — s) (s) ds] dv. (17) 


To determine the condition for the existence of a finite, mean value in the stationary condition(t—> ) 
we divide the interval of integration with respect to y into two sub-intervals, Let us pick some value of T 
which is finite but sufficiently large for TR,(T) «0. Then we can set approximately 


( R,(s)ds = (s) ds, sR, (s) ds = (s) ds. 
0 0 0 0 


We can write 


m (t) [— mv + o3|@—s) R; (s) ds| dv + 


mal ox {—(m,— alr, (s) (s) ds} do. 


Frein here it follows that the mean value for the stationary condition will be finite if Condition (16) is 
satisfied, 


Letting the coefficient of mutual correlation Ry2(r) = 0 in Equation (14), we obtain the equation for the 
two-dimensional moment 


My (ty, ts) — v) + exp (ts + tz—u—v)+ 
0 ; 
+508 (ty, a) +408 (ta, + 


+ (ty, (ta, R, (%_ — da, dz,} du dv. (18) 
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Without giving the coefficients of correlation Ry and R, concrete values it is difficult to find the conditions 
for which, for example, the dispersion in the stationary condition has a finite value, 


It should be noted that particular cases of the example in question were analyzed in papers [2, 3]. 


In paper [2] the one-dimensional probability density for the process y(t) was derived for zero initial 
conditions, with y(t) determined by the differential equation dy/dt + € (t) y = € (t)h, where hg is a unit step 
function, 


The results of paper [3] are applicable to the example under consideration if m,=0, Let us show this. 


It follows from Expression (17) that in this case m(t) =0 and Equation (18) gives the correlation function 
which can be written in another form, Let us introduce instead of the correlation function of the fluctuations 
n(t) their spectral intensity S_ (w) by means of the well known formula 


(u — v) = dey, (19) 


—oo 


Let H(jw, tedwt be the response of System (1) to the harmonic signal ce". For zero initial conditions 
we have from Expression (2) 


H (jo, t)= exp (jou mt dz) du. 


It is not difficult to verify that, after substituting Expression (19), the exponential multiplier appearing 
under the integral sign in Equation (18) can be represented in the form of the mean value of the product 
H (jw, ty) H(— jw, ty). Then Equation (18) can be written in the form 


(tas ta) Sa(w) (Je, ty) H (— ta) do, 
For ty =t, tg=t+T, we obtain the equation 
k(t, = (0) <H (jo, t) H (— jo, t + do, 


which is given in paper[3], From here it can be seen that the last equation is of a purely formal nature and 
does not reveal the calculation difficulties, 


In concluding we will note that the formulas obtained above can be generalized for the case when the 
correlated, random functions €(t) and (t) are not normal but possess a probability density in the form of many 
. dimensional Edgeworth series [4], However, the final formulas and calculations turn out, in this case, to be 
very cumbersome and unwieldy and are. therefore not given here. 
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“A METHOD OF SOLVING A PARTICULAR CLASS OF INTEGRAL EQUATIONS 
BY CASING COMPUTERS 


Iu. S, Validenberg 


(Moscow) 


Zeidel's iterative method of the approximate solution of Friedholm's and Volterra's 
first and second kind of linear integral equations is presented, The feasibility of using 
computers for the solution of special types of these equations is indicated. The applica- 
tion of this method to some problems of synthesis of automatic control systems is examined, 


INTRODUCTION 


In the solution of problems of synthesis of automatic control systems (ACS), to find their optimum 
characteristics it is necessary to solve integral equations [1], The necessity of mechanizing the solution of 
these equations stimulated the development of a method for their solution using computers. The essential 
nature of such a project is also determined by the fact that a large number of problems in many fields of 
science and technology leads to integral equations, In the present work a method is given for the solution, by 
means of a computer, of Friedholm's and Volterra's first and second kind of linear integral equations, having 
a convolution-type kernel 


b x 

9 (2) = f(z) K (x—8) p(s) ds, (1) 
b Or x 

| K(e—s)e(s)ds, (2 


where the given functions are the kernel K (x— s) and f (x) the right side of the equation, g(x) is the solution, 
0 =x,s=b, The given functions are relationships obtained experimentally and introduced into the machine 
without a preliminary approximation, The problem of finding a unique solution is solved. 


The integral equations are solved approximately by the method of successive approximations, i.e., the 
iterative method [2, 4]. 


Method of Solving Equations 


The iterative method of solving integral equations is used not only in numerical solutions but also in 


special computers, This method is very convenient for use in computers because of the uniformity of the 
mathematical algorism. 


Two basic types of the iterative process are known: ordinary iteration and Zeidel's iteration (or Gauss- 
Zeidel) (2, 3], In the case of Zeidel's iteration, when the m-th approximation for the i-th ordinate of the sought 
after function is calculated, the earlier found 1-st, 2-nd,... , (i 1)th ordinates of the m-th approximation 
for the sought after function are taken into account; while the i-th, (i+ 1)-th, etc, ordinates are taken from the 
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preceding (m— 1}th approximation, In other words, refinements of the results of the calculations are carried 
out after each ordinate of the sought after function is determined, and not after the calculation of the given 
approximation for the whole sought after function, as is done in ordinary iteration, This explains why Zeidel’s 
methods leads, in many cases, to a much more rapid convergence than the one arising from the ordinary method 
of iteration, 


Up until now, the method of ordinary iteration [4, 6, 8, 9] has been used in computers developed for solving 
integral equations, Zeidel's iterative method has been given preference in this paper, Besides speeding up the 
convergence, as mentioned above, this type of iteration permits one to construct an iterative process in a form 
which secures the convergence of Friedholm's first kind of equation with a symmetrical kernel, the quadratic 
form of which is positively determined, In the case of ordinary iteration convergence of the process may not be 
guaranteed for these equations, Besides this, Zeidel’s iteration, as will be shown later, leads to simpler computer 
circuits than those arising in ordinary iteration. 


The given integral equation is approximated by a system consisting of 32 linear algebraic equations, This 
system is solved by Zeidel's iterative process, The suitability of the approximation for equations of the 2-nd 
kind is obvious, while for equations of the first kind it is valid in the case of symmetrical kernels, the quadratic 
form of which is positively determined, It is easy to show that for such kernels the solution of the integral 
equation is unique (if it exists), which permits one to approximate the integral by replacing it by a finite summa- 
tion, 


The approximating system is solved by the iterative method using the following algorism [3]: 


. n—1 
= fi—0,5 AsKy ol? — As — As Kuo, (3) 
(4) 

in? jm 

f ee b 

= Kuri (ims, 2,.3, As => i n=32). (6) 


Here the index in the raised brackets indicates the number of the iterative step, Kj; is the matrix approxi- 
mating the kernel K(x— s), Equation (6) is associated with Volterra’s equation of the first kind, and represents 
a recursive relationship realizing the method of successive substitutions, 


Functions arising in integral equations are, as a rule, continuous and have continuous derivatives, In 
isolated cases where this is not so, special artificial methods have to be used. 


A zero vector is used as the initial approximation of the iterative process, In those cases where the zero 
vector leads to a diverging process (although these cases are possible, they are relatively rare) the integral 
equation is solved in two stages [4], First, a rough and very approximate solution of the integral equation is 
found by the minimizing method (i.e., along an open cycle) in accordance with the following algorisms: 


2 
fi —0,5 AsKi,9, — As » Kis; = ¢?-» min, (7) 
n a-1 
py [ — As » K jp; — (AsK + 1) | = min. (8) 
im] 


In the second stage a more accurate solution of the integral equation is obtained by the iterative method, 
with the rough solution, obtained in the first step, being used as the initial approximation, The introduction of 


an extra step greatly increases the method's possibilities since it permits one to solve, for example, Friedholm's 
equation of the second kind for cases where a parameter of the equation is greater than the minimum eigenvalue, 
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As an illustration of the nature of the convergence of the iterative process we give three numerical 
examples of equations solved in accordance with(3) - (5). 


0,9 
(2) = e-* (x + 1.5) —0.08265 — | e 9(s) ds, 


0 
e(z)=e*, n=30, k=6, 


9 (2) = — 9 (2) = 5.24e n= 30, 


=0, 0.6, k =8, 
= 1 +\e*9(s)ds, 9(z) =0.5 (1 + e%), n= 30, = 0, k=5. 
0 


Here k is the number of iterative steps, The error in the solution of these equations does not exceed 1%, 


Experimental Verification of the Convergence of the Iterative Process 


One can find in the literature descriptions of several devices for the solution of integral equations realizing 
the method of simple iteration [4, 6, 8, 9). These computers contain, as a rule, two basic circuits: the cal- 
culating unit and the memory unit (see, for example, [9]). The calculating unit carries out che multiplication, 
integration and subtraction operations in accordance with Equations (1) and (2), The memory unit memorizes 
the ordinates of the sought after function which are determined in the calculator, and which are thereupon re- 
introduced into the calculating unit, In this way, the iterative process is realized with a closed system con- 
sisting of two units ~ calculating and memory. 


In the case of ordinary iteration the memory device consists of two memory units, one of which stores the 
k-th, while the other stores the (k + 1th approximation of the sought after function, After the calculation of 
each approximation, a special switching device switches these units in such a way that the (k + 2)-th approxi- 
mation is recorded in the unit in which the k-th approximation was stored [9]. 


One, and not two, memory units are required for the realization of Zeidel's iterative method and, more- 
over, the switching circuit is not required. 


Any of the known methods of analog computation can be used in the calculating unit of the device. Thus, 
for example, standard electrical mode! units were used to this end in [9]. 


Operational memory devices, found in general purpose digital computers, can be used as the memory unit. 
The memory cores in such devices usually have two stable conditions, which guarantees the stable operation of 
the whole device, The location device can be replaced by a simpler arrangement which guarantees the successive 
scanning of the cores in the memory device. 


The following two problems [5] can be considered to be the basic problems which arise in the develop- 
ment of a computer for the solution of integral equations. 


1. The problem of designing a memory device, Because Zeidel's iterative method is used, he problem 
can be solved with relatively simple circuits, 


2, The problem of designing the generator of a kerne! which is a function of two independent variables. 
In the general case, for an arbitrary form of kernel, the problem is soluble but leads to very complicated cir- 
cuitry. For convolution-type kernels it is not necessary to design special equipment, and one can make use of 
a conventional generator of a function of one argument; while the convolution is achieved by the proper dis- 
tribution in time of pulses which control the operation of the apparatus, In this case, experimental relationships 
are used as the given functions in the integral equation, The fact that the problem is soluble only for convolu- 
tion-type kernels has been completely justified in practice, since many problems of analysis and synthesis of 
ACS which reduce to integral equations, lead to a convolution-type kernel, 
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The principles which form the basis of the method (convergence of Zeidel's iterative process, stability of 
the closed system, etc.) were verified experimentally on a special apparatus, the schematic of which is shown 
in Fig. 1. 
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Fig. 1. Circuit of a unit for the verification of the convergence of the iterative 
process: 1) motor, 2) trigger circuits, 3) integrator, 4) storage device, 
5) kipp oscillator, 


Here the persistence of a cathode-ray tube is used as the memory unit, The function being recorded is 
fed to the deflection plates of the tube. The trace of the beam does not disappear at once and the image of the 
function is read by a photomultiplier FEU-1 with the help of a Nipkov scanning disc with holes distributed in a 
spiral, 


Photomultipliers FEU- 2 and FEU-3 are used to obtain the horizontal sync and vertical sync pulses, The 
integral transformation computer (ITC) integrates the function under the integral sign. The control pulse circuit 
(CPC) consists of frequency dividers and a kipp oscillator, and produces the pulses which ensure the normal 
operation of the equipment, 


x 
The apparatus was used to solve Volterra’s equation of the 2-nd kind (xz) = 1 —|j @ (x — s)ds. 
The solution g(x) =e"* was obtained in 2.5 sec with an error of 5-6%. 


Possible Applications of the Computer 


In the synthesis of ACS it is necessary to find the characteristics of a system satisfying the given technical 
conditions, The transfer function @ (s) or the pulse transient function k(t) is such a characteristic, In the 
solution of this synthesis problem the structure of the system being synthesized is usually not known beforehand 
and, therefore, the corresponding differential equation cannot be set up. 


The characteristic of a linear system which realizes the transformation of definite actions, is given by 
Volterra’s equation of the first kind with a convolution-type kernel 
t 
z(t) =\ g(t — (e)de, (9) 
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The synthesis problem in the indicated statement is very actual for the case when the system is acted on 
by signals which are random functions of time. In this case the synthesis problem includes the problem of find- 
ing the optimum (in a definite sense) characteristics of the system. Thus, for example, if the acting signals are 
stationary random functions, given by the corresponding correlation functions R(r), the k(t) of the system, 
resulting in the minimum mean square error, is given by Friedholm's equation of the first kind [1] 


| Re (A) dh = Rag (2). il 


Equations (9) and (10) can be solved on the equipment, with experimental relationships used as initial 
conditions, As an example, let us examine the special, but very interesting from the practical point of view, 
case when white noise is taken as the interference, In this case the filtering problem is reduced to Friedholm's 
equation of the second kind 


k(t) = (t) — | Rm (t — dh, Ry (*) = y°8(*), 


The solution of this equation for R,,(r) = e"® J I c? = 8 (see the second example, see p, 714) is obtained 
in eight iterative steps with an error of 1%, Equation (11) can also be solved automatically, without the 
participation of an operator; at the same time the voltage, corresponding to R,,(r), acts from the correlator, 
while the quantity y? is expressed as a level of the signal's 
spectral density, measured in that part of the spectrum 
where the useful signal is small, 


Methods of solving Equation (10) for other concrete 


sin wt cases may require additional refinement. 
janet Fields of possible application of the equipment are 


not limited to the solution of integral equations with a 
conyolution-type kernel and with experimental curves 
(e.g. correlograms), It is also possible to solve integral 

equations with analytically defined kernels, These kernels 
are realized on analog computer amplifiers which permit repetitive operations at a frequency of 50 cycles, 
while the control pulses, determining the work cycle, are formed in the control pulse circuit (CPC). Figure 2 
gives an example of the use of such amplifiers for the calculation of Fourier transforms, The kernels cos wt 
and sin wt are obtained as solutions of the determining differential equation x" + kx = 0 [7]. This can be used, 
for example, when it is required to ‘ind the transfer function of a system and the solution of the integral equation 
gives the pulse transient function, ‘fhe Fourier transform of k(t) permits one to obtain the transfer function 
directly, 


SUMMARY 
1, An apparatus is developed for the solution of a definite class of integral equations by means of Zeidel's 
iterative process. 


2, Provision is made for the solution of integral equations with experimental curves and convolution- type 
kernels, as well as with analytically established kernels realized on an analog computer amplifier. 


“In conclusion I must express my thanks to V. V. Solodovnikov for his guidance during the execution of this 
project. 
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AUTOMATIC OPTIMALIZER 


A. A. Fel'dbaum 
(Moscow) 
The paper deals with the problem of constructing a machine which minimizes a 
function of several variables Q (x4, . . . , Xn) in the presence of the additional restrictions 


Hj(x4,...,%,) <0 (j=1,...,m). Several possible ways of solving this problem are 
examined, and the circuits of different models of the machine are presented. 


1. Statement of the Problem 


Self-adjusting systems, which are at present in an early stage of development, will apparently become 
very important in the next few years. 


One of the most important classes of seli- adjusting systems are automatic optimization systems, the pur- 
pose of which is to find the conditions for a minimum in some quantity Q, usually under additional restrictions. 
Let us examine the case when Q is a function of the variables xy, ... , Xp, and firstly the variables x; (Fig. 1, 


4) are set up at the input to the unit O, after which, following a certain elapsed time, the quantity Q appears at 
the output of the unit: 


Q=Q (ti, stn). (1) 


The function Q is not known beforehand and the values of x;, which minimize Q, are unknown. 


In Fig. 1 the apparatus A automatically selects the quantities xy, .. . , X, so that a minimum in the 
quantity Q would be obtained for the general case— with the existence of the additional restrictions, in the 


form of inequalities, which the variables xj must 


a 
> > A; <0 
a Aa a 
: The number of restrictions m is arbitrary; the 
;|< i < kf three possible cases are m < n, m=n, and m> n, An 
apparatus which carries out such an operation is sub- 
sequently called an automatic optimalizer, In the 
Fig. 1 absence of restrictions the output Q of the unit O is 
connected to the input of the automatic optimalizer A 
Fig. 1, a). If restrictions do exist, the boundary values Hj(j=1,...,m) are also fed from the output of O 


to the inputs of the automatic optimalizer A (Fig. 1, b). 


Let us assume that Q has only one minimum or that the initial conditions do not originate in a definite 
region of attraction to the one minimum of Q, In the presence of several minima the described apparatus be- 
comes more complicated; however, it is not difficult to extend the principles being examined to the case of 
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several minima, Also, let the functions Q and Hj be continuous and differentiable over the whole region of 
variations of x;.* 


Automatic optimalizers can be utilized as machines for the direct solution of variation problems, for the 
automatic synthesis of optimum control parts of automatic systems, for the automatic determination of the 
dynamic characteristics of complicated (for example, nonlinear) members, as well as for the optimization of the 
behavior of controlled objects under production conditions. 


2. Possible Methods of Designing an Automatic Optimalizer 


The basic method of solving the given problem is the determination of a method of automatic hunting. 
In[1, 2] some of the hunting methods were examined: the method of successive changing of the variables, the 
relaxation method, etc. The following methods are effective for the machine described below. 


a. The gradient method, in which the gradient vector is first determined: 


n 
gradQ@ = >) > ij. (3) 
j=1 
Here i; are the unit vectors in the space x4, ..., Xp, Where the x; are assumed to be rectangular co- 
ordinates, After the gradient is found at the given point, a step is made in the direction opposite to the direction 


of the gradient, i.e., in the direction of the quickest decrease in the magnitude of Q, Following this the gradient 
is again determined at the new point, another step is taken, etc. 


b, The method of steepest descent (3j in which, after the gradient is determined, motion occurs in the 
direction opposite to that of the gradient until the minimum of Q is found along this direction, Then, an analogous 
cycle is carried out, consisting of the determination of the gradient and the subsequent operational motion, etc. 


After considering the possibilities of the different methods with respect to the rate of convergence, as well 
as the convenience of introducing various restrictions Hj, for the case of appreciable distances from the minimum, 
preference was given to the method of steepest descent, This method has the following advantages: a) speed of 
convergence, especially significant when the number of variables x; is large; b) realizable by means of a re- 
latively simple circuit; c) convenience of imposing any type of restriction, 


The last condition is very important. For example, in the application of the method of successive changing 
of variables (the Gauss- Zeidel method) the circuit is simple only if there are no restrictions, Their presence 
greatly complicates the circuit, 


The imposition of restrictions is convenient if the gradient method is used (see [4)), but this method takes 
second place to the steepest descent method on the basis of the speed of convergence, Because of this it is better 
to use the steepest descent method where large deviations from the minimum are involved, However, for small 
deviations, the gradient method yields the minimum, stable, mean deviation from the minimum and guarantees 
better tracking of the slowly shifting minimum. Because of this, when small deviations from the minimum are 
involved, operations in the apparatus being described are based on the gradient method, Circuits for the method 
of steepest descent and the gradient method are similar, and transition from one method to the other is no problem. 
We can find analogies to such a combination of methods in some automatic control systems where different laws 
of motion are selected for large and small deviations so as to obtain a combination of short regulation time in 
the transient process, and small error in the steady state. 


3. Versions of the Optimalizer without Restrictions 


As an illustration of the machine's operation, when large distances from the minimum are involved, let 
us examine the simplest case of Q as a function of two variable x, and xq (see Fig, 2 where the (x, x3) plane is 
shown), Let the initial values of xy and xg correspond to the point My on this plane, Each cycle consists of a 


* Incidentally, this condition is not necessary. It is sufficient that discontinuities of the first kind do not exceed 
some limiting value. 
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*trial* motion followed by a ‘working’ motion, The trial motion consists of a series of consecutive steps Ax; 
in the direction of each of the variables separately, and the determination of the components (9Q/ xj) of the 
gradient grad Q at the point My. The working motion is in the direction opposite to that of the gradient, At 

the same time the variables x; change by discrete increments 


gradQ cus Ax=—a (4) 

& J [23 where a =const > 0. In this way, a transition occurs from 

ft Pig. point Mp, along the straight line MyM, to point My where 

Sf the function takes on its minimum value along the given 

Fg A * straight line [5], At this point another trial movement occurs, 

$4 j M the gradient is determined, and this is followed by a new 
/ iz working motion along the straight line MyMg, etc. This con- 
ian tinues until the system arrives close to the minimum and, for 

gradQ 


example, the quantity = Q/ax| turns out to be 
Fig. 2 


smaller than some preassigned, small constant quantity £4. 
If £3 <€ <€, (where €, is a fixed quantity) then the automatic search can be carried out in small steps with 
small values of a (accurate search), If € > £ , then the search is carried out in large steps (coarse search), 
Various sizes of steps make it possible to obtain various combinations of accuracy and speed of the automatic 
search, Each step, whether trial or work, is accompanied by a measurement, and occurs over a time interval 
called the access time. 


Determination of 
gradient 


| 


ccurate ste 
land nd check of fc 


step 


one 


Fig. 3 
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There are two possible versions of the system: 


a, alternative coarse and precise scans are made along each line MjMj,4,: having located the minimum 
approximately, the system uses small steps to locate it more precisely; 


b. a constant size of step is used, determined from the value of € before the commencement of each 
work motion, 


The principle of operation of version "a" is clarified with the help of a block diagram of the sequence of 
operations, shown in Fig. 3, The operation of determining the gradient, with which the machine's work cycle 
starts, is shown at the top of the block diagram, This is followed by a selection of one of the next possible 
directions of the operation (operation E), If, for example, — > €, then the next operation is a coarse step; if 
E; <€ =< €z, then an accurate step is made; finally, if € = €4, then no working motion ensues and, after some 
time 1, the operation of determining the gradient is repeated. 


Each coarse step is followed in turn by a selection— operation E, depending on the sign of the increment 
Rev — of the magnitude of the i-th step in question, where 


AQ = — Qu—n- (5) 


Here Q;) is the magnitude of Q after the i-th cycle, while Qi; -, ) is the same quantity after the (i— 1)-th 
cycle, If AQ < 0 then the coarse step is repeated; if AQ > 0 then a reversal occurs and subsequent steps are 
accurate ones in the direction opposite to the initial one. Operation E is repeated after each accurate step. If 
AQ <0 then the accurate step is repeated; if AQ > 0 then a reversal and one accurate step in the new direction 
(which now coincides with the initial direction 
of the coarse hunt) are made, With this the 


f cycle for locating the minimum ends, After a 
Determination of time r a determination of the gradient is re- 
gradient peated, 


The circuit may be simplified by omitting 
F ri A the last two steps in the cycle enclosed by a 

: broken line; however, the accuracy of the search 
| is reduced in this case, 


bbe L_5<5 Figure 4 is a block diagram showing the 
E sequence of operations in version “b", It is 
simpler than the preceding one since here there 
| } 7 { is no reversal, The size of step on each line 
coarse step and accurate step a MjMj 41 depends on the value of which emerges 
icheck of AQ check of AQ from the operation of determining the gradient. 


Otherwise the designations and nature of the 
operations are here the same as in Fig. 3, and 

E E therefore do not require special explanation, 

The search accuracy in this system, for a given 

4 number of cycles, is smaller than in the preceding 


a<o| [aaa] [agro 
L__ The automatic optimalizer consists of two 
parts— the operational and control units, The 
Fig. 4 first can be realized using analog-type calculating 


devices (abbreviated CD), the second using digital- 
type CD, Thus, the automatic optimalizer appears as a combination of a digital and analog computer, However, 
no digital-to- analog, or vice versa, converters are necessary since the analog section represents a closed loop 
which is the operational unit, The digital section control unit— receives at its input binary quantities (i.e., 
quantities having only two possible values) from the analog section and, in turn, controls the analog section by 
means of switches. The principle of operation remains unchanged, irrespective of what type of equipment is 
used — electromechanical, electronic, etc, As an example, Fig. 5 shows the schematic of a completely electronic 
operational unit for the “a version of equipment. 
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Fig. 5. 
Before the unit is put into operation, the zeros are set up on the dc amplifiers Ag - Ag,..., Asn ~ Asn» 


Aq and Ag, and the initial conditions are set up on the integrating sections— amplifiers Ag; - Agn. The voltages 
Ugy > Ugn at the output of these sections represent, on a particular scale and with a reverse sign, the initial 
values of x;, which are delivered to the Sitput of unit 0 through the inverters Ag - Agn. After the zeros and 
initial conditions are set up, the start button (on the control unit) is pushed and the circuit begins to operate 
cycle after cycle, A large cycle begins with the operation of determining the gradient, This operation in it- 
self consists of a series of access cycles, 


1, Memorizing the “old” value of Q, For this, the electronic switch Ky at the output of the unit O, and 
switches K, in the circuit of amplifier Ay, which are normally closed, open for a time r and then reclose.* 


* The opening and closing of the switches is done by the control unit. 


122 


A A 
| 
8 
I 
| Ale Ul V ie 


At the output of the unit O there is already, at this time, a value of Q corresponding to the values of x; delivered 
to the input, If switches Ky are open, which can be represented by short circuits,* then section Ay is inertial, 
Consequently, its output voltage will become, after a known time r», equal to the quantity Q with a negative 
sign, i.e., equal to Qy = — Q(x, ... » Xp). 


2, Memorizing the partial derivative 8Q/8x,. Here switch Kg, through which a small increment Ax 
is transmitted to the input x, of unit O, is closed. Unit O receives a ‘start’ signal from the control unit. The 
unit goes through a work cycle (for example, a transient process occurs in the electronic model) as a result of 
which there appears at the output of this unit, after some time, the value Q, = Q(x, + Ax, Xg,..., Xp). Now 


switches K, are closed but switches Ko, Ky, Kaj are opened for a time r, As a result, Ag develops at its output 
the voltage 


= Qo = Q + Ln) + = — AQ). (6) 


When switches Ky and Kj are open, section Ag, represents an inertial section with an amplification 
factor equal to unity, Therefore, the voltage Ugy = AQ, is established at its output and is subsequently 
memorized there when Kg and Ky are closed, In view of the constant nature and small value of the incre- 
ment Ax, the voltage Ug, is approximately proportional to @Q /dx,. 


3. Memorization of the partial derivative 9Q/axg. This is done in an analogous manner, Here switch 


Kj} is opened, the increment Ax is delivered to input x,, as a result of which Ag develops at its output the 
voltage 


Us = Qo — = La, Ln) — Q(t, + Ax, In) = — AQ,. (1) 


Then switches Ko, Kya, Kgg are opened for a time r and the quantity AQs, proportional to 9Q/ dx, is 
memorized in Ag». 


Analogously, the memory circuits Ag - As, are caused to memorize the values AQ, ~- AQn, which are 
proportional to the derivatives 8Q/ ax, - 8Q/dx,. With this the operation of determining the gradient ends, 


After the end of this operation a selection is made, from the value of €, of the nature of the subsequent 
operations (Fig. 3). To obtain the value of € , with various weighting coefficients,** for |@Q/ ax; |, 4 well 
known circuit containing diodes can be used, Let us assume, for example, that > E,, Then a coarse cycle, 


consisting of a series of small cycles, goes into operation, Each small cycle consists of the following access 
cycles, 


1. The ‘old’ value Q; is memorized, At the same time all the switches Kg - Ksn open for a certain 
time, Then the quantities 9Q/ax; multiplied by certain coefficients, proceed to the inputs of the integrating 
sections Ag; - Asp. Consequently, the quantities xj; receive incremental changes determined by Formula (4). 
One of the groups of switches, Kg - Kgn or Kj - Kh, can be eliminated from the circuit, The operation of 
one or the other group depends on the sign of the Ax delivered to the input of O, 


2. A “start” signal is given to the unit O, and a work cycle occurs, As a result there appears, after some 
time, at the output of unit O a new value =Q + AXxgj, + AXei, . » + AXn 


3, The control unit analyzes the sign of the difference AQi+ =Qi+, — Qi, which is formed at the output 
of amplifier Ay, with switch Ky open, All the switches in the memory circuit Ag - Asn are closed at this time. 
If AQi4, <0, then the first of the above cycles is repeated. If AQi4, > 0, then (Fig. 3) a reversal occurs and 
subsequent steps become smaller in size, i.e., they become “accurate”, 


In the circuit this is achieved very simply. Instead of switches Kg, - Kg, opening simultaneously with 
every step, switches Kj - Kj}, now open simultaneously and deliver to the integrating sections voltages opposite 
in sign to the former ones. If the new input voltages are delivered through larger resistances (Ry > Ry) then the 
changing currents of the capacitors in the integrating loops will be smaller, Hence, the voltage increments at 


the outputs of the integrating loops Agy - Agn, also decrease by a corresponding number of times, i,e., steps Ax, 


* In the laboratory model, series-type, diode, two-pole switches, controlled by triggers, were used, 
** If weighting coefficients need be taken into account. 
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will become smaller, Resistors Ry and R, can be made equal but, to obtain a large stop, two switches Ky - Kjy, 
for example, would be opened simultaneously, while for a small step only one of these keys would be opened. 
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The same cycles are repeated for an accurate search that are used in a coarse search, When AQi+, 
becomes greater than zero, then a reverse and one accurate step backwards must follow. This is achieved, for 
example, by switching on, for one cycle, switches Kj - Kj,, which deliver, through larger resistors Rg, to the 
inputs of the integrating sections, the voltages Uy; - Ugn which are opposite in sign to the corresponding voltages 
Ugy ~ respectively. 


At the end of the operations described above the gradient is determined again and a new large cycle 
begins. 


let Ex £4. In this case the control unit does not begin a new large cycle, but, after some time, repeats 
the operation of determining the gradient. 


The described principle of operation makes it possible to obtain not only two grades of search steps, but 
also any number of such grades, for example ‘coarse’, ‘medium’, and ‘accurate’ steps, The principle of 
achieving such a circuit is apparent from what has been said, It is also possible to control the step gradually 
depending on the value of €. However, experiments have shown that there is no real need for this. — 


Figure 6 shows a simpler schematic of the basic unit, corresponding to the "b” version and the block 
diagram in Fig, 4. Since there are no reversals here, invertors Ag - Agn, which were present in the first version 
(Fig. 5), can be omitted, 


As in the block diagram in Fig, 4, the large cycle here begins with the determination of the gradient, As 
a result of this operation, voltages Ug - Us, become proportional to the partial derivative 9Q/ax,- 8Q/dXp, 
respectively, 


The analysis of the quantity € is carried out in the same way as in the preceding circuit, If this is 
followed by an accurate search, then switches Kg, - Ksy are simultaneously switched on for a certain length of 
time, if it is a coarse search then switches Ky - K',) are also switched on, Let us note that in one channel it 
is sufficient to have only one switch, say Kgj, but it must be turned on for a longer time if a coarse search is 


_ tequired, and a shorter time if an accurate search is required, In this case the circuit of the operational unit 


is somewhat simplified, but the control unit circuit becomes more complicated, 


The analysis of the sign of AQ; can be carried out, for example, by means of a simple circuit, the output 
of which can have only two values: upper (nominally 1) and lower (nominally 0), When U,= — AQ; > 0, the 
output of the circuit is in state 0. When U,= — AQ < 0, the output changes to state 1. 


Fig. 7 


Switches Ky - K'yn can be omitted from the circuit in Fig, 6, and one grade of search step is sufficient 
because the automatic optimalizer includes an automatic step control, In fact, it can be seen from Equation 
(4) that during the work motion, the increment Ax; is proportional to | @Q/ax;j| . For appreciable distances . 
from the minimum, i.e., where | 0Q/ax; |is usually large, the size of the step Ax; is also large. Close to 
the minimum, where | 8Q/8x; | is small, the size of the step also decreases, Because of this the laboratory 
model of the apparatus did not contain a changeover from a “coarse” to an “accurate” search although, in some 
cases, such a transition is extremely desirable, The quantity — was used for another purpose — for the change- 
over, for small €, from minimization by the method of steepest descent to minimization by the gradient method, 


Figure 7 shows the block schematic of one version of the control unit, For the sake of simplicity, the 
control unit for the simplified circuit in Fig. 6, for n = 2, is examined, However, analogous principles can be 
applied to the design of a control circuit for a more complex system, shown in Fig, 5. 


Triggers T;, Tg, Ts and Ty, are connected in a closed loop, When trigger T, is transferred from state 1 
to state 0, its output pulse passes through the delay circuit r and through either switch Kp or Kjgy , and trans- 
fers trigger T, from state 0 to state 1. Trigger T3, in its turn, returning from state 1 to state 0, switches over 
trigger Ts; from 0 to 1. In the same way, trigger Ts acts, through switch Kogr, on Tg, and T, acts on Ty. Besides 
the large loop Ty, Tz, Ts and T,, a small closed loop can be formed from T, and T,, if switches K, and Kim 
are open, 


Trigger Ty gives the ‘start’ signal to unit O, when it changes from 0 to 1, and the ‘stop’ signal, when it 
returns from 1 to 0, Trigger T,, operates (changes to 1) for large € (E = € ), while trigger Tog, operates for 
small values of € <€)). <€ triggers T,, and Tog are in the 0 state, 


Trigger T;,,, operates when AQ> 0 and is in state 0 when AQ < 0, 


Let us trace the sequence of operations in the control unit, Switch Ky is in the output of the pulse generator 
PG, which supplies low frequency pulses, Switch Ky is only open when all the triggers Ty - Tg are in the 0 state. 
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In this case (as well as when a trigger pulse appears) the pulse arrives at the 1 input of trigger T, and transfers 

it into state 1. As trigger T, operates, it opens switches Ky, and K', (Fig. 6) as a result of which a value of Q 

is memorized. After some time, a pulse from generator PG;, operating at an appreciably higher frequency 

than PG, (these generators do not have to be synchronized with each other), appears at the 0 input of trigger Ty. 

Consequently, T, switches over to state 0 and its output pulse passes through either of the switches Kp and Kim: 

and operates trigger T,, This trigger opens switches Ky, K', and Kj, Here, the quantity x,, with increment 

Ax is fed to the input of unit 0, the quantity Q changes and the increment AQ, proportional to 9Q/ax, is 

memorized by the section containing amplifier Ay. Then trigger T, releases under the action of the next pulse, 

which arrives at its 0 input from generator PG;, and causes trigger Ts to operate, This trigger opens gates Kg», 
‘mn and Ky, and the quantity 9@Q/8x, is memorized in the circuit of amplifier Agg, Then trigger Tj, is 

put into the 0 state and trigger T, operates — the ‘work’ motion trigger. It opens switches Kg, - Kgg (Fig. 6), 

because of which x, and x, undergo incremental changes by Ax, and Ax, in accordance with Equation (4), After 

one step, trigger T, releases and trigger T, operates, 


The subsequent sequence of events depends on the magnitudes of € and AQ. If € = £2, then trigger Tp 
is in the 1 state, switch K. is closed, and Rp is open, Then trigger T,, as it releases, causes trigger T, to 
operate, Another work motion occurs, These motions continue until the quantity AQ becomes greater than 
zero, Then T;,, operates, switch Kir closes, and Kj;, opens, Now the operation of determining the 
gradient begins, i.e,, the sequence of operation of triggers T,;, T, and Ts, If € = € 2, then this operation is 
followed by a new work motion, i.e., the sequence of operations of Ty, Ty, etc. If, however, € <€,, then 
trigger T,., Operates, switch K,,, closes and opens the loop of triggers. They all return to state 0, but sub- 
sequently, a pulse from PG, excites the loop after some time through switch Ky. 


If &; < & < €y then the search operates according to the gradient method and not the method of steepest 
descent, 


In fact, in this case trigger T is in state 0, switch K, is closed, and K, is open, This is followed by 
the consecutive excitation of triggers Ty, Tg, T, and Ty, i.e., the determination of the gradient, and one work 
step, This is followed by a new determination of the gradient and another work step, etc. Trigger T, operates 
in the case when at least one of the triggers T;, T; and T, switches over to the 1 state. Then the ‘start’ signal 
is given to the unit 0 which can be, for example, in a particular case, an electronic unit, The next pulse from 
generator PG, returns trigger T, to the 0 state ("stop” signal) before triggers Tz, Ty and Tg switch over to the 
0 state. Triggers Ty, T:, Ts, Tg and Ty can be replaced by monostable flip-flop oscillators, This simplifies 
the circuit, since generator PG, is no longer required — the circuits return themselves to the 0 state, However, 
this is accompanied by complications in the resetting of the unit's operating speed. 


4. Automatic Optimalizer with Restrictions 


The principle of operation of this optimalizer is illustrated in Fig. 8 for the case of two variables x, and 
Xq and two restrictions H, < 0 and H, <0. However, this principle is applicable with no change to the case of 
any number n’ variables and an arbitrary number m re- 


strictions, 

H-0 N The lines H, = 0 and H, =0 are plotted on the (x,, x») 

“grad @ plane, and the ‘forbidden’ regions H,> 0 and Hg> 0 

‘ie 40 \ NN are shaded, While the objective point (xq, x,) moves in 

PMT NK same way as the one discussed before. However, 

YW ma yO N even one of the quantities Hj becomes greater than zero, 
YY the system's operation changes. The machine's ‘attention’ 


Bs ee goes the same operations which previously the quantity Q 
as underwent in the permitted zone, Suppose, as an example, 
the quantities Hjy, Hjz,..., Hjz (1 < m) become greater 
than zero, Let us denote 


A, = Ay +-Higt+...+ 
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Now the operation of determining the gradient proceeds as follows: using the same method as before, 
dH, /ax; (x; =1,..., 0) are calculated and are memorized at the outputs of the same sections Ag - Asn 
(Figs, 5 and 6) as before, Then the objective point moves along a line ( — grad Hg) until all the quantities 
Hy, (k=1,...,2) become negative, As soon as this is achieved, the unit returns to analyzing the quantity 
Q alone, 


In Fig. 8 the objective point moves along the line MpMy, i.e., along the line (— grad Q), Im the process 
it enters the forbidden region H, > 0 and arrives at point Mj, Hence, further on, Hy is determined and the point 
proceeds along the line M4Mg, i,e., along the line (— grad Hy), At the point Mz the magnitude of H, becomes 
less than zero, Consequently, further motion proceeds along the line (— grad Q), etc, 


Figure 8 shows, on an exaggerated scale, the motion of the objective point along the boundary H =0, 


Various cases of a minimum on the boundary are possible: either the function Q has a minimum at the 
point Ny, on the boundary Hy = 0, or the value of Q decreases up to the point Ng— the points of intersection of 
the boundaries H; = 0 and Hg=0, In the first case, the objective point, upon reaching N,, begins to perform 
small oscillations about this point, or else stops completely if € is small, In the second case the point crosses 
the line Hy =0, in which case Hy and Hg become positive and the system begins to act on the quantity H, = H, + 
+ Hg, and leaves the forbidden zone along the line (— grad Hy), Here the objective point begins to perform 
small oscillations in the vicinity of the point Nz. With this the search ends since Ng is the required point, At 
the point corresponding to the minimum, located on the gradient boundary, Q and Hy are directed toward each 
other, 


For an optimalizer with restrictions it is necessary to add to the circuit described earlier an additional 
device, the circuit of which is shown in Fig. 9. 
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When any of the quantities Hj, at the output of the unit 0, become positive, the triggers corresponding 
to them, from the set Ty, = THm? operate, Dc amplifiers can be introduced in front of these triggers to in- 
crease the system's contin to the sign of Hj, When trigger T}4; operates, it opens switch Ky, at the input 
of amplifier A";, As a result, invertor A} develops at its output the necessary sum Hg, in accordance with 
Equation (8), The operated triggers also act on the diode circuit D, - D,, (an "Or" = type logic circuit) at the 
output of which appears the positive voltage Uj,, The voltage Uj, acts through the switch K,,) and the 
cathode- follower on the trigger TQH: which opens the switch Ky; and closes switch Ky, Now Ho, and not Q, 
is passed to the input of amplifier Aq (Fig. 5). Beside this, a negative pulse, from any anode of any of the triggers 
Ty,, proceeds through the diodes Dj, Dj - Dp, Djp to the input of the trigger T,,., (Fig. 7) and puts the control 
unit into the “determination of gradient" condition, Let us assume that the system begins to operate with the 
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quantity Hy. As soon as the operation of determining the gradient begins and trigger T, operates (Fig. 7), it 
causes Tp) to operate (Fig, 9). The latter closes switch Kj, and the state of trigger Tc}; “freezes* until the 
end of the operation of determining the gradient. At the end of this operation trigger Ts, in releasing, switches 
trigger T),, to the 0 state, and switch Kp) opens. As soon as all the quantities Hj become negative, all the 
triggers TH - THm release (return to the 0 state), the quantity Un, becomes zero and, after Kj, is opened, 
releases trigger Toj;. Now the quantity Q is once again passed to the input of amplifier Ay. The pulse from 


the triggers Tjj, causes trigger T,., to operate, In this way, the circuit is once again returned to the condition 
of “determination of gradient." 


There are other methods of including the restrictions, with which the search rate along the boundary is 


increased somewhat, However, this is accompanied by appreciable complications in the circuit of the restrictions 
unit (Fig. 9). 


For large values of n it is convenient to introduce certain changes in the optimalizer, permitting thereby 


a great economy in components, The problems of designing such systems, for large values of n, must appear as 
the subject of a separate paper. 


The systems considered above are not universal, In practice there arises the possible case where, instead 
of an isolated maximum in the space (x, x2, . . . , Xp) there is an (n— k)-dimensional space, where k =n— 1, 
+++» 1, in which Q has the same maximum value [6]. In this case the logic of the system's operations to 
achieve such a sub-space must be more complicated; here a much larger number of trial movements becomes 
necessary, However, experiments have shown that even the relatively simple logic described above, leads to 
success in many important examples, The problems of complicating the logic for the special cases mentioned 
above, as well as to allow for the effect of variations in the function Q(x, ... , Xp) itself during the time of 
the trial motion, if this does occur, are the subject of a separate paper. 
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TWIN-CHANNEL AUTOMATIC OPTIMALIZER 


R, I, Stakhovskii 


(Moscow) 


The circuit and some units of a twin channel, electronic, automatic optimalizer 
for the location of minima (taking restrictions into account) are described in detail, 
Experimental results, corresponding to its different operating modes, are given, 


INTRODUCTION 


The object of this work, some of the results of which are presented here, was the design of a machine for 
the location of the extremum of a function of many variables, based on the main systems suggested in[1], The 
design of the optimalizer required both theoretical [2] and experimental study, This paper is devoted to the de- 
scription of a twin- channel optimalizer which permitted us to solve the following problems, 


1. To verify, experimentally, the principle of operation of hunting systems, suggested in [1]. 


2, To conduct an experimental evaluation of various hunting algorithms and to select the best operating 
algorithm for the apparatus, 


3. To develop different types of simple and reliable circuits for the different units in the machine and to 
plan ways of further simplifying them and increasing their reliability, 


4. To verify experimentally the theoretical relationships between the hunting accuracy and time and the 
system's parameters, for different methods of hunting. 


5. To determine experimentally the optimalizer's characteristics for its operation with a standard single 
capacity plant characterized by the equation Q = Ayx} + A,x3 Here a study was made of plants with constant 
and variable parameters and with added restrictions, A study was also made of the condition under which the 
plant was acted on by random interferences, 


Systems for the automatic location of extremums, with or without added restrictions, can be used as the 1 
basis of automatic machines designed for a wide variety of applications, beginning with a machine for the 
solution of variational problems, and ending, for example, with an automatic-operator in a production plant. 


The design of a twin-channel type of optimalizer represents only the first step in the development of 
multichannel systems, Consequently, the scheme was designed with the possibility of it being developed for 
the case of a function of n variables in mind, This resulted in the choice of a circuit with a small number of 
amplifiers per channel, which makes the scheme adaptable to an increase in the number of variables, To 
simplify the circuitry it was necessary to select a simple operating algorithm for the optimalizer, from the types 
described in [1], and a simple block diagram of the control unit, in which the hunting method can be changed 
during the actual hunting. 


Some of the problems listed above are clarified in[1), Therefore, there is no need to dwell on them 
here, Let us note only that the methods of simplifying the over-all system, as well as its separate units, are far 
from exhausted, Further simplifications must be made during the design of multichannel systems. 
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General Description of Equipment 


The complete schematic of an optimalizer which will hunt when two restrictions are imposed, is shown 
in Fig. 1. The equipment is divided into three units: operational unit, which is the nucleus of the equipment, 
control unit and restrictions unit, The details of the block schematic are described in({1) and, therefore, only 
a general description of the equipment, as well as separate, particular additions and changes, are given below, 


Fig. 2 


Since a simplified algorism, in which there is no reversal, is used in this equipment, each channel con- 
sists of three amplifiers with the necessary number of electronic switches, The input of the operational unit is 
a device which determines differences between the values of the input function at two discrete instants of time 
("differentiator"), This device has two inputs: for the minimizing and resulting restriction functions H», The 
setting-up of initial conditions is carried out on the step integrators A, and Ag, Trial increments are given by 
means of switches in the summing amplifiers Az, Ag. 


The control unit is designed as a two-circuit loop, composed of the single flip-flop oscillators Oy (O,H), 
Oz, Os, Oy, which control the group of switches Ky(Kyy), Kz, Kg, Ky, in the operational unit, 


The switching of the loop circuit, which is necessary when a change is made in the hunting method, as 
well as when the steepest descent method of hunting is used, is achieved by means of switches Kp, Rp. Kim 
Kim» Which are controlled by triggers Tp and Tim. 


To determine the instant at which the hunting algorism is switched, a unit is used which measures the 


2 
quantity |8Q/ Ax, |. 

The restrictions unit is designed to switch the inputs from Q to H», by means of the transfer of control 
from switches Kyo to analogous switches Kyy, if even one of the quantities Hj (x;, xg), i=1, 2, would exceed 
zero, At the same time, the operation of determining the gradient is initiated. In order to avoid false steps, 
trigger To};, which switches the inputs, is locked during the determination of the gradient, The unit contains 
triggers TH and Ty, with amplifiers Ay and Ag, The latter are provided for more precise operation of the 


triggers. Tu, and Ty, are designed to introduce one or two restrictions by means of switches Ky, and Ky,, 


as well as to connect in the device for determining the gradient, every time the boundary is crossed. The pulse 
generator PG serves to ensure the uninterrupted operation of the loop of monostable flip-flops, 


Figure 2 shows a photograph of a laboratory model of the optimalizer. The operational unit is in the center, 
The operational amplifiers are situated in the lower part of the front panel, Above them is the memory circuit 
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with lugs to which the inputs and outputs of the amplifiers are brought out, The input resistors and those used in 
the feedback loops are situated in the memory circuit or else in the switching circuits on the rear panel of the 
unit, 


The test panel, with meters xq, xp, Q(x;, x9), AQ, £, is situated on top, In front of the test panel, and 
set up on the operational unit, is a model of a nonlinear unit with a quadratic characteristic, Next to it, on 
the right side, are manual controls for setting up the initial conditions, 


The control unit is situated to the right of the operational unit, and is equipped with a memory circuit, 
This circuit is used to display the operation of the switches by means of neon bulbs with light when the switches 
are ‘1 a conducting state, Various parts of the control unit, in the form of circuit sub- assemblies, are distributed 
on the front and rear panels of the unit, 


The restriction unit, located to the left of the operational unit, is also equipped with a memory circuit 
for its own operations, Amplifiers Ag, Ayg and Ay, Ayg, triggers TH» TH,» TqH and the corresponding electronic 


switches, are situated on the front and rear panels of the unit, 
It should be noted that experimentation has shown that the use of display-type memory circuits is very 
expedient, because it helps in the adjustment and checking of the equipment's operation in different modes. 


Short Description of Separate Components 


1, Operational amplifier. The requirements for an operational amplifier used in the optimalizer, 
differ appreciably from the requirements for computing amplifiers used in computers. Actually the require- 
ments for accuracy and permissible amount of zero drift are less stringent, but the reliability requirements are 
higher, especially if industrial applications are contemplated, Actually, the amplifiers in the optimalizer are 
servo-amplifiers, and errors in them lead merely to a slowing down of the hunting process, but do not affect the 
value of the determined minimum, The only exceptions are the amplifiers in the input section of the operational 
unit, which demand more precise adjustment, However, even here the circuit could be simplified considerably. 
A very simple and reliable circuit of a three-stage, operational amplifier, using 6N9 and 1/2 6N8 tubes, and 
deriving its power not from a vibrapack but from regulated power supplies (+ 300, — 190, — 300 vy), is utilized 
in the apparatus. Zero drift over an extended period of time (weeks), does not exceed + 1.5 v, which turns out 
to be quite sufficient for the reliable operation of the whole equipment if the range of output voltages is + 100 v. 


2. Electronic switches, Two-pole, series switches are used in the operational amplifiers, The circuit of 

such a switch is shown in Fig. 3. When a positive voltage, with respect to ground, is delivered to point 1, si- 
multaneously with a negative voltage to point 2, diodes 

D, D, and D, conduct since positive voltages, with respect 
to cathodes, appear on the anodes of both diodes, The 
output signal passes through the switch by way of D, or 
D3, depending on the polarity of the signal, When the 
controlling voltages on points 1 and 2 are reversed, 
diodes D, and D, disconnect the input resistors from the 
input to the amplifier becuase the anodes of D, and D, 
become negative with respect to their cathodes, due to 
the low forward resistance of Ds and D,, 


The switch's most important characteristic, as it 

is used in the optimalizer, is the small drift of the 
operational amplifier, with electronic switches in its 
input, when the amplifier operates as an integrator. The 
drift is caused by a parasitic emf which exists in the 
equivalent circuit of the switch, and which depends, all 
other things being equal, on the cathode temperature 
and the dissimilarity of the parameters of diodes Dy and 
D,. In a circuit containing an integrating section, which 
operates from the switch, it is easy to obtain a drift not 
k= kilchm, m = megohm exceeding (when the section's time constant Rin C= 1 sec) 

Fig. 4 + (1 - 3) v/min for the open and closed states of the switch. 
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3, Controlling monostable flip-flop, The electronic switches are controlled by monostable flip-flops 
in the control unit, which are circuits with one stable state, The delay time is controlled by the grid bias, 
The circuit has a push-pull output, A ring of flip-flops is obtained by connecting the input of a 
multivibrator to the corresponding output of the preceding one through a capacitor of 50 -100 uf. whereby the 
return of one monostable multivibrator to its stable state would trigger the following multivibrator. 


4, Finder of € as a sum of the moduli of partial derivatives, The circuit for determining € is shown in 
Fig. 4. The arrangement operates in the following manner, 


If the input voltages are positive then the gates, shown on the left in the schematic, conduct positive 
currents through the resistors R and the summing resistor r, If the input voltages are negative then the gates on 
the right form a sum of negative voltages on the second resistance r, The difference between the voltages on 
the resistors r, proportional to the sum of the moduli of the input voltages, is extracted by means of a subtracting 
stage and is amplified by the required amount in an amplifying stage with a cathode- follower output, 


5. Model of a nonlinear unit with restrictions, To test the optimalizer, a nonlinear converter, for two 
input quantities, was designed as an analog to a unit characterized by one minimum and approximated by the 
equation 


Q = Ajai + Agr}. 


Part of the circuit of the nonlinear element, shown at the top of Fig, 5, is similar to the circuit of the 
Em finder if the input resistors are replaced by thyrite resistors, having approximately a square law characteristic, 
The subtraction stage is replaced by two operational amplifiers, To minimize the loading of the last stage of 
the amplifier the thyrite resistors are followed by current dividers, 


The nonlinear element's characteristic is shown in Fig. 6, During experimentation, the resistance of the 
summing amplifier's feedback loop was shunted by a capacitor to simulate the unit's model's time constant 
Tout ¥ 9.7 sec, It should be pointed out that the optimalizer being examined was intended primarily for 
operation with purely discrete units (units of the 2-nd kind), and the experiments were conducted with a model 
of a unit of the 1-st kind [3] only for the sake of simplicity. 


The introduction of restrictions was simulated by means of the circuit shown in the lower part of Fig. 5. 
Here the boundaries of the regions H, (x;, xg) > 0 and H,(x;, x,) > 0 are straight lines, having a positive and a 
negative slope with respect to the abscissa in the plane (xy, xg), which can be moved in this plane by changing 
the values of voltages Ay and Ag, 


Results of Testing of Automatic Optimalizer 


Object with constant parameters. Figure 7 shows the dependence of Q(t) on time when the steepest 
descent method of search is used. The search begins with a determination of the gradient which is marked by 
two overshoots of AQ on the horizontal plateau Q (x49, Xy). Then steps are taken in a direction opposite to that 
of the gradient until the next value of AQ becomes positive. After this the gradient is again determined (two 
overshoots on the horizontal plateau) and more steps are taken, etc. until the minimum is located, 


It can be seen from Fig, 7 that despite the relatively fast first stage of the search (already in the first 
step Qmin © Qpre» insofar as in the object Ay w Ag, i.e., the gradient at any point in the plane (xq, xq), is in 
the direction of the minimum) the total search time is nevertheless large, because an appreciable length of 
time is required to deplete the error which is obtained when the minimum is overshot due to the large size of 


step. 


Figure 8 shows the diagram of a hunt based on the gradient method, Here overshoots occur at each step, 
and the value Q.,,,. is approached very smoothly, The stabilized value of the mean excess fg over Qpre (fy = 
=Qmean in the stable state) in the case of the gradient method of search is considerably smaller than for the 
steepest descent method, 


Figure 9 shows the search time T», and the mean excess in the stable state f, over the magnitude Qpre» 
as a functions of the scale of step for a steepest descent method of search, The scale of step was varied by 
beams of the amplification factor of A, in the input differentiator. From Fig. 9 it can be seen that the curve 
Ty (k) has a minimum. This can be explained by the fact that for a small scale of step the search time is large 
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because of the slow motion of the objective point, while with large scale steps a lot of time is lost in depleting 
the errors in determining the minimum, For large scales of step, the mean excess fy is practically independent 


of k, For smaller scales of step, fy increases because of the operating threshold present in the trigger for detecting 
the minimum 
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Figure 10 shows the dependence of the search time T, and the mean excess fy on the size of step scale for 
the gradient method of search, where the scale of step was varied by means of the trial step x, 
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The broken lines are the experimental curves Fig. 10 
calculated from the formulas 2): 


In (z, 43) 
infi fo = 0-2582*(A, + Ay), 


where Xo,,,, is the initial value of the maximum coordinate (xomax = 60 v), 6 is the magnitude of the zone 
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around the minimum, in which the transient process of the search is considered to be finished (§ =6 v), At is 
the time taken to make one step (At = 10.5 sec), a is the proportionality factor between the size of the incre- 
ments in the coordinates and the magnitude of the partial derivative (a = 2), Ay =» A, =A is the coefficient of 
the squared term in the object's equation, A = 0.0215, §x is the size of the trial step. 


The figure shows good agreement of the theoretical and experimental curves Ty (k), if one takes into 
account the fact that the object's characteristic is not an exact quadratic, 


Comparison of the search diagrams shows that the steepest descent method leads to a more rapid approach 
to the minimum than does the gradient method, However, close to the minimum, the search time with the 
chosen method of locating the minimum is practically the same as the search time with the gradient method. 
At the same time, the accuracy in determining the value of Q,,,, in the case where the gradient method is used, 
is considerably greater than in the case of the steepest descent method, because of the existence of an operating 
threshold in the minimum- locating circuits. 


_ Thus, it is expedient, in setting up the first initial search step, to utilize the steepest descent method, 
while, when in the close vicinity of the minimum, which can be judged from the value of €, the gradient method 
should be applied, 


Object with variable parameters. Experiments with an object with variable parameters were conducted 
with the gradient method of searching, since, as was already pointed out, it is necessary to use this method of 
searching in the vicinity of the minimum. 
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Fig. 11 


A very important property of the optimalizer is its ability to track the moving extremum point in the case 
where its coordinates are functions of time. For example, let this function have the form: 


Q = A, + (t)}? + Az [%2 + 
where @4(t) and @,(t) are the given functions of time. 
It can be shown that there exists a critical value 


a of the derivatives *,(t) and 6‘,(t) when the optimalizer 
céases to track the minimum effectively, and the search- 
5 ing becomes disorganized, since the system begins to make 
false steps. 
In fact, between the time when the function 
Xg9) is memorized and the instant when the quantity AQ 
, / b is memorized, there elapses a time At. If, during this 


time, #4 (t) and $,(t) change sufficiently, then the 


f <9 the object's parameters are independent of time. A critical 


speed will be attained when the magnitudes and signs of 


v the derivatives and ',(t) are such that AQogr, 
a 4005 bbtfrey when decreasing to zero, changes sign, This is obtained 
Fig, 12 


effective value of the increment AQgfr, which is memorized 
in terms of the values of the partial derivatives, can be larger 
or smaller than the true value of AQ, which is obtained when 
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when 


©; (t) = 82/At (i =1, 2), 


where §x is the size of the trial step, At is the time interval between the memorizing of the quantities Q (yp, 
Xm) and AQ). 


As a result of this situation, when the coordinates of the minimum oscillate with a frequency Q( 4, (t), 
@,(t) are sinusoidal functions of time with amplitudes Ax, and Ax,) the output, when in a stable state, also 
contains harmonic components with frequency Q, 20, 3Q, etc. A plot of the output value of Q, obtained when 
(t) =11 sin and (0) = t, is shown in Fig, 11. 

Figure 12 shows a plot of f, as a function of frequency of oscillation, for constant amplitude of oscillation, 
for the case of oscillations along one coordinate only (curve a), and for the case of simultaneous oscillations 
along both coordinates (curve b) with §6x =15 v. The mean deviation was determined as the ordinate of a 
rectangle with its area equivalent to the output oscillations, 
In case b, a more rapid growth of f, occurs; in both cases, 

f the rate of growth of f, increases as the critical speed is 
. approached during a certain part of the period, 


Ae, =11258 One way of increasing the critical speed is to increase 
i. Ae, 8.75 8 the trial step §x, Figure 13 shows a plot of f, as a function 
of the size of step §x when the object swings along both co- 
rad ordinates, and the scale of step is kept constant, From the 
ag: oh graph, it is evident that when §x is decreased f, increases 
rt rapidly as the critical speed is approached, The relationship 
a $ 0 Sz.y between fy, in a stable state, and the size of the step’s scale 
Fig, 13 is interesting in the case of sinusoidal oscillations of the 
minimum, with §x=const, The graph in Fig, 14 shows that 
in this case, there is an optimum scale of the step, for which 
fo» V the quantity fy is a minimum, In fact, when the scale is 
j small, the system is unable to follow the oscillating minimum, 
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while with a very large step the system overshoots the mini- 
] mum. 


&r=65 ¥, / An object with random interference, A series of ex- 
\ / periments to determine the realizability of a search in the 


\ presence of random interference were carried out to determine 
6 7 the stability of the optimalizer in the presence of interference, 


The interference was created by means of a noise generator, 
developed by A, M, Petroyskii [4], which produces approxi- 
Y . , 4 ok mately “white noise over the range of frequencies of interest 
Fig, 14 to us, namely, a few cycles and lower, The noise, together 
with the optimalizer signal, were fed to the input of the ob- 

ject, so that the signal at the object's input appeared as the sum of the trial step §x and the random interference, 

the scale of which could be altered by changing the generator‘s output voltage, In the experiments the noise 

signal was varied from zero to rms values exceeding, by about a factor of two, the value of 6x, Comparisons 

ipa made by recording on the potentiometer EPP-09 the composite signal's coordinates x; + 6X + fran. noise 

=1, 2). 


Figure 15 shows the dependence of f, and Ty, on the output voltage of the noise generator, As can be seen 
from the graph, it can be assumed, as a rough approximation, that f, is proportional to the noise level. At the 
same time, To is relatively independent of the noise level, This situation is explained by the appreciable 
suppression of interference due to the persistence of the memory loops, These, being high-pass filters, neutralize 
noise signals having a zero constant component, but are very sensitive to the quantity §x, i,e., to the constant 
component of the composite signal, For objects discrete in time this type of interference suppression is impossible 
because the interference comes through in the variation of the quantity AQ, which appears as a constant quantity 
at discrete instants of time, When operating with discrete objects, it is necessary torepeat the supplied §x and to 
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average the resulting series of values of AQ. Figure 16 shows a graph of the search x, = f (t) in the absence of 
noise and in the presence of very strong noises, 


Object with added restrictions, In a hunt, in the 
presence of added restrictions, two situations can arise: 
a) the minimum is in a “permitted” region, b) the 
minimum is on the boundary of a forbidden region, 


In case “a” the search proceeds in the same way 
as in the absence of restrictions, except, maybe, in the 
case where part of the path is along the boundary of 
the region, The location of a minimum situated on 
the boundary is the more interesting case, It was in- 
dicated in[1) that, in this case, the stabilized process 
oscillates in the vicinity of that point where the gradients 
of Q and Hy are directed exactly toward each other, 
Figure 17 is a diagram of the transient process when the 

minimum in the permitted region is on the boundary, 
2U, 


Figure 18 shows the plot of a search for the case when 
the minimum occurs at the point of intersection of two 
Fig, 15 boundaries, If in the first case (Fig. 17) the amplitude 
of the steady-state oscillations is constant, then in the 
second case (Fig. 18) overshoots of various sizes occur. 
They occur when the operating point is pushed out of the region where both restrictions are violated simul- 
taneously and the absolute value of the gradient Hy is approximately doubled. 


It is apparent from the diagrams that the descent along a boundary during the transient process is, natural- 
ly, slower than an unrestricted search, the reduction in search speed being by approximately 1.5-2.5 times. 


SUMMARY 


1, The described instrument, designed on the basis of the principle proposed in [1], can be used to locate 
the extremum point of a function of several variables, the shape of which may not be known. The requirements 
imposed on the object the extremum of which is being located are merely that this extremum — actual or con- 
ditional (restrictions taken into account) — exists. 


search in lsearch in absence’ 
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Level of minimum without restrictions 


Fig, 17 


2, The best algorithm for such a device is a combination of the method of steepest descent and the gradient 
method, It has been shown experimentally that it is possible, in the determination of Qmin(Cpre)» to attain an 
accuracy of a few tenths of one percent with respect to the maximum possible value of Q, 


3. The developed circuits of the elements of the equipment turned out to be satisfactory from the point of 
view of simplicity and reliability, but there is still room for further refinements both in the over-all structure of the 
apparatus and in the circuits of the separate units, 


4, Experimental verification of theoretical relationships (where this was possible) showed satisfactory agree- 
ment between the experimental results and those obtained-from design formulag derived in [2]. 


5, The experimental study of the equipment's characteristics showed that the parameters of the optimalizer 
which have the greatest effect on the hunting process are the amplification factor of the input differentiator and 
the value of the trial step § x, as well as their product, proportional to the scale of the work step. 


Whether the object's parameters are variable or constant, there is a definite, optimum scale of step, which 
gives the minimum value of the mean deviation in the stable state, For each value of §x there is a maximum 
permissible speed of motion of the extremum point, The direction of the critical velocity depends on the direct- 
ion of the tial step § x, and the permissible volocity increases with increase in §x, It should be pointed out, how- 
ever, that a small change in the method of determining the gradient, namely the allowance for the time derivative 
Q'(t), permits a considerable increase in the permissible velocity of the minimum's variations, 


A study of stability in the presence of interference showed that, when objects which. are continuous in time 
are involved, the device copes with interference quite successfully, Here the signal-to-noise ratio reached 1 :2 
at least, An increase in the noise level results in an increase in the search time and the amplitude of the stable 
oscillations, but does not interfere with the operation of the equipment, An increase in the stability in the 
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presence of interference can be achieved by a repeated feed (not one at a time) of a trial step with a subsequent 
averaging of AQ. 


The author thanks the director of this work A. A. Feldbaum, and A, V. Kalinina,who participated in the 
e perimental part of the work, 
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NEW PRINCIPLES OF CONSTRUCTING TELEMETRY SYSTEMS WITH 
PULSE-TIME AND PULSE- WIDTH MODULATION 


V. A. Il'in and A, I, Novikov 


(Moscow) 


Considered are the principles of constructing single-channel and multi-channel 
telemetry systems with exponential transducers which allow creation of multi-channel 
systems with time separation of channels without the use of commutators, 


The basic expressions for these transducers are given, as well as examples of 
systems, 


One of the principal special features in the telemetered transmission of information is the linear depend- 
ence of the measurement erro: on the stability of the transmission coefficients of the entire system and of its 
individual elements, A variation of the transmission coefficient by § % gives rise to the same error. Therefore, 
the problem arises in telemetry of creating high-stability transducers (modulators, demodulators, etc.). 


B 


Fig. 1. Circuit of an exponential pulse-width 
transducer, 


As a result of work done on this problem at the 
Institute of Automation and Remote Control of the 
AN SSSR, new, simple, high-stability pulse-time and 
pulse-width transducers, called exponential, were 
suggested and developed, This paper describes such 
transducers and their underlying theory, and considers 
examples of telemetering systems. 


The exponential pulse-width transducer con- 
sists of a bridge (Fig. 1), in one arm of which are connected resistors Ry and Rg, and in the other arm of which 
are connected a resistor and a capacitor or a resistor and an inductance, A diode forms the bridge's diagonal, 
In the circuit of Fig. 1, Ry « R and the back resistance 6f the diode, Rqp > R, If a rectangular voltage (Fig. 
2a) is impressed on the bridge, the voltage on condenser C varies as shown in Fig. 2b. 


Fig. 2, Time diagrams of the voltage and current 
in the transducer, 
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At the moment of time t,, the voltage on the condenser reaches the voltage on resistor Rg and the diode 
begins to pass current, Thus, from the moment t, until the end of the pulse, the voltage on the condenser re- 
mains constant, since the direct resistance of the diode, Rgg « R, R > Ry 


Figure 2c shows the variation of current through the condenser or, alternatively, the voltage on load 
resistor R;. 


a 


‘Fig. 3. Circuits for exponential transducers. 


The width t, of the pulse at the output load does not depend, with this circuit, on variations in the source 
voltage: 


te = RC in 


Choosing as our data unit either resistor R or capacitor C, we obtain a linear exponential transducer, in 
which the dependence between the pulse width t,. and the data unit resistor R or capacitor C will be linear.* 


In many cases it is expedient to use as a data unit a low-ohmic resistor Rg (Fig. 3a) or to control the 
transducer by feeding the input voltage to the diode (Fig. 3b). With this there will hold a logarithmic de- 
pendence between the pulse duration and resistor Rg or the transducer's input voltage, However, by using a 
similar circuit to carry out the reverse functional transformation at the receiver, one can obtain a linear de- 
pendence between the magnitude of the original parameter and the parameter at the receiver output. Such an 
exponential transducer will be functional, An example of such a system is given in Fig. 6. 


If the internal impedance of the pulse generator is taken into account, the pulse width. t, will be 


t=C(R+KR,) 
| 
where 
~ Ri +R," 


For « andR » +R,, 


| 


If the internal impedance varies by ARj, the relative error is 


— 
RF ER, 


* The idea of a linear exponential transducer was suggested by V, A, Il"in. 
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The absolute change in pulse width for Rj « R and R » Ry+R, will be 


The characteristics of actual diodes differ from the idealized piece- wise linear curves, so that the diodes 
begin conducting current earlier than the moment of time t,. This leads to a broadening of the trailing edge 


of the pulses and to a dependence of the pulse width on the operating threshold Ae of the diode and the shaping 
circuit (Fig. 2b), 


Ra, [le 
% 
Pa, 
fa, 
R 
Fig. 4, Transmitter cincuit for a three-channel Fig. 5, Time diagram of the currents and voltages in 
telemetering system. a three-channel transmitter system. 


The absolute change in pulse length when the threshold varies by Ae, will be 


Os, = —C(R+ KR,)In 
Ri + 


With actual elements, a transducer can be constructed with a repetition period T from several tens of 
microseconds up to tens of seconds with an error not exceeding 1% and a trailing edge of less than 1%, * 


Pulse-time modulation is implemented by introducing a differentiating network at the transducer output 
and pulse-shaping stages in the transmitter, 


Besides the well-known methods of implementing multi-channel systems with time and frequency separa- 
tion, there is a good deal of interest in methods of constructing multi-channel systems with pulse-time modula- 
tion and time separation of the channels without the use of a commutator, Figure 4 shows a simplified trans- 
mitter circuit for such a three-channel system, Here, Ra,, Ra, and Ra, are potentiometer data units of the first, 
second and third channels and resistors Rg, ,Rqg and Rg, are used for creating the shielding time intervals 
between the channels, The differentiating network R},Cj is connected to the circuit output, 


* The idea of the reverse functional transducer was suggested by A, I, Novikov. 
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We note that if, at the input of the transducer channels, a voltage is applied, then the circuit is set up in 
accordance with Fig. 3b. 


Figure 5 is illustrative of the maximum time deviations in the channels and of the shielding intervals 
between channels, 


1 

< 
2 2 
¥ 


Fig. 6, Circuit for a telemetering system with exponential transducers. 


Figure 6 shows a circuit for a transmitter and receiver (on the right) of a single-channel VST-1 tele- 
metering system, developed at the Institute of Automation and Remote Control of the AN SSSR as a first attempt 
to replace the antiquated Bristol Company system and its modifications. Here, 1 is an indicator, 2 is an 
electronic relay and 3 is an audio oscillator. The transmitter system is constructed from semi- conducting 
elements and from other elements of high reliability (Fig. 7), It contains no vacuum (radio) tubes, contact 
elements or local power supplies. A photograph of the transmitter is shown in Fig. 8. 


Fig. 7. Transmitter circuit of a telemetering system. 


The de pulses generated in the receiver have a width of 1.5 seconds and a repetition period of 2 seconds 
(Fig, 2a), They are fed to the transmitter system by line. The working pulse voltage at the transmitter input | 
varies from 25 volts to 60 volts, with an input impedance of about 1200 ohms, The system was designed for 
standard potentiometric data units with impedances from 300 ohms to 2000 ohms, 


_ The electronic relay, composed of germanium triodes, is connected to the output of the exponential 
transducer, The output triode of the electronic relay is simultaneously an audio oscillator, with a constant fre- 
quency chosen within the limits of 200 cycles to 2000 cycles. 


At the instant t, (Fig. 2), the electronic relay operates and switches in the audio oscillator until the end 
of the initial pulse, 
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At the receiver input is a resonant circuit, tuned to the frequency of the audio oscillator, The voltage 
in the resonant circuit is amplified, and causes the relay to operate at the instant t,, The relay switches con- 
denser C into the output stage for a brief period, and the output 
stage then “remembers” the condenser voltage until the following 
cycle, 


The RC networks in the transmitter and receiver have 
identical time constants, This guarantees a linear relationship 
between the magnitude of the measured quantity and the indica- 
tion of the output instrument (an inverse functional transforma- 
tion). 


The distinguishing characteristic feature of the receiver is 
the inclusion of a memory element, If the state of the data unit 
changes, the indication of the output instrument remains constant 
until the next measurement cycle (until the moment t,), after 
which the next indication is established during the time constant 
of the instrument, 


Tests revealed the following characteristics of the system, The measurement error corresponds to a class I 
instrument. The additional errors arising as the ambient temperature of the transmitter varies from— 50° to 


+ 50° and the supply voltage varies by + 15% are less than 1%, The radius of action, with the equivalent of 
four- millimeter bi- metallic lines, was greater than 500 kilometers, 


The system works normally when the series resistance in the line varies from 0 to 10,000 ohms, The 
frequency stability of the audio oscillator was + 1.5% when the temperature and supply voltage varied within 
the limits stated, 


It should be said, in conclusion, that exponential transducers, and telemetering systems with exponential 
transducers, should find widespread application in various branches of science and technology. 


Received April 1, 1958 
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NOISE STABILITY OF TELEMETRY SIGNAL TRANSMISSION 
ON CHANNELS WITH NOISE FLUCTUATIONS 


Vv. A. Kashirin and G, A. Shastova 


(Moscow) 


The potential noise stability of transmission parameters is determined under ampli- 
tude, frequency, pulse-frequency, pulse-time, pulse-width and pulse-code modulation. 
Generalizations of the transmission parameters are introduced, It is shown that in some 
range of the signal-to-noise ratio, the very best noise stability is attained with frequency 
modulation, 


INTRODUCTION 


We consider a transmission parameter varying within the limits + 1 in a random way, but such that, 
with a sufficiently high accuracy, it can be represented by its values at moments of time which are separated 
by a fixed time interval T, The interval T is related to the highest frequency Fy, of parameter variation 
by the well-known relationship 


T a) 


We shall assume that all parameter values in the interval + 1 are equally probable. 


Let the parameter be transmitted along the 
channel by means of a signal A(,t) which, in the 
general case, is a function both of time and of the 
transmitted parameter. The function A(, t) depends 
ne on the type of modulation, Figure 1 shows the time 

variation of the parameter \ and the function A(A, t) 


for pulse-width modulation, 


Fig. 1 channel with weak noise fluctuations, the specific in- 
€ tensity of which equals g volts/cycle, errors arise, the 
relative value (over the entire range) of which may be determined, for an ideal receiver, from the formula 


_ 


#0.) = (2) 


Formula (2) was obtained from the corresponding expression in[1] for the absolute value of the mean 
square error by dividing the error by the magnitude of the range of variation of the transmitted parameter and 

by changing the limits of integration, The integration of the square of the derivative @A(A ,t)/dA is carried 
out in the present work over the limits 0 to T, while in [1] the integration was over the segment from — T/2 to 


T/ 2, Such a change of limits is completely admissible and is not essential from the point of view of the theory 
of potential noise stability, although it does influence the resulting values of the mean square errors for frequency, 
time and width modulation, 


The mean square error over all values of A is found from the usual integration 
+1 
\ 8 (Ag) (3) 


-1 

By means of the formulas given, we now determine the potential noise stability of various forms of modu- 
lation, with the condition that the dynamic range of the signal AQ,t) is limited, i.e,, the signal can vary 
within the limits — Uy, to Um or 0 to 2U,. 
Amplitude Modulation (AM) 


With amplitude modulation, the signal can be written in the form 


t) =U, (4 md) cos apf, 


where U, is the mean value of the signal amplitude, m is the index of amplitude modulation and wp» is the 
signal carrier frequency, 


For finding the error, we define the integral: 


Iau = dt= cost gta. 
0 0 


Iam = for T 
By substituting the value of I, 4 in (2), we obtain the formula for the error with AM: 
bam 4U m*T : 


For convenience in comparing the various transmission methods, we introduce the generalized parameter 
p, equal to the ratio of the maximum signal voltage to the effective noise voltage in the frequency band taken 
by the parameter A: 


i+m 


By taking into account that with AM, U,,, =U, (1+.m), we get Say = 


The minimum error is for m =1, 


Thus, for amplitude modulation, the error is completely determined by the parameter p and increases in 
inverse proportion to it. Formula (5) remains valid for any noise level [1]. 
Frequency Modulation (FM) 


Using frequency modulation, with carrier frequency fy and deviation frequency fg, the signal can be pre- 
sented in the form 


A(h, t) = Umsin + 


8am min= V 2/0. (5) 
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T 
Um (28)? \e cos* (fy + tdt. 


If T then Ipy = /, Ut and the error is 


We now introduce the parameter y, equal to the ratio of one-half the band taken by the signal to the 
frequency band of parameter A: 


(6) 


= (7) 


By taking (1), (4) and (7) into account, we obtain the expression for the mean square error using frequency 
modulation: 
t= @) 


In contradistinction to AM, the error using FM depends not only on parameter p, but also on the quantity 
Y, characterizing the degree of band expansion needed by the signal in comparison with AM. As is obvious 
from (8), for y >» 1, the error, with weak noise, increases in inverse proportion to the first power of the broaden- 
ing coefficient of the frequency band. 


The use of secondary AM impairs the noise stability [2]. 
With this, the mean square error is 


2.2 


Pulse-Frequency Modulation (PFM) 
A, Pulse- frequency modulation with constant pulse width (PFMD. 


With PFMI, the frequency of the sequential pulses f;Q.) which have a constant form and amplitude, 
varies linearly as a function of the magnitude of the measured parameter, \: 


AO) 


where f; is the mean frequency of the pulses of the sequence and fg is the deviation frequency. 
The signal can be written in the form of a sum 


where F(t) is a function describing the envelopes of the individual pulses and ty is the moment of transmission 
of the k’th pulse, | 


k 
ETN (k= 1,2,...). 


| The integral sought is 
| 0.146 
Af 
With frequency modulation, 44 =fqg+F,, and 
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The derivative @A(A ,0/ dA will equal 


8A (a, t) OF (t —t,) 
In order to compute the derivatives 9F(t— t;)/9A, we introduce a change of variables: 


Since the pulses are nonoverlapping, and the derivative dF(t;,)/ dt, # 0 only during the time when the 
k'th pulse exists (where the magnitude of the derivative is independent of the number k), the desired integral 
equals 


We now consider transmission by means of triangular pulses of amplitude 2U,, and width r, 
In this case 


z= 


By the use of (9), (10) and (2), we find the magnitude of the squared error: 


30 (fy + 140 
128 U3, T + 1.5L + 0.5) 


== 
where L = T(f; + fd) is the number of pulses appearing during time T when the given parameter A is trans- 
mitted, 


As is obvious from the expression obtained, the magnitude of the squared error depends on the values 
assumed by the parameter, By remembering that all values of A within the limits from — 1 to+1 are equi- 
probable, we determine, from Formula (3), the mean square relative error, averaged over all values of A: 


30% 1.5 


+ In + — in [0.5 + 1.5 7y + — 


1.42 
where x =f; + fg and y fy. 
fee 1 (or f; > fq), i.¢., if during time T for any A, a large number of pulses is transmitted and r= 


a , then the expression for the squared error simplifies to: 
8 


PEMI 428 U2, ae 


Then, 
(ium 1602, 
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Taking into account that Afppyy = 2f4y, and making use of (4) and (6), we obtain an expression for the error 
in generalized parameters: 


(13) 


On the basis of (11) it may easily be shown that the error is significantly decreased if the deviation fre- 
quency is increased, With given frequency limits f, and f,,, the maximum value of f, is attained for 


(4) 


he 

It will be assumed in the sequel that Condition (14) always holds. In the case when the frequencies fm 
and f, are not given, but only their difference Af =f,,— f,, is given, the problem arises of choosing the working 
portion of the frequency spectrum, By using (12), one can assure oneself that the derivative 06 /df, is 
positive, Consequently, from the point of view of noise stability, it is advantageous to decrease f,, given the 
band Af, If the highest frequency of variation of the transmitted parameter is F,,,, then the least value of f, 
equals 

Formula (12) is valid, however, only for fy, f; > fy. 


To investigate the dependence of § on f, for small fy, we return to general Formula (11), By taking 
fy = k(2F,,,) and fg >> fj, we find that even for the low-frequency portion of the spectrum, 06 / dk > 0, i.e., 
the error also falls when the lowest frequency of the range is decreased. 


Thus, the minimum error occurs when 
fi = 2F m (15) 


By substituting (15) in (11), and using generalized parameters, we obtain an expression for the minimum 
error (in the case when f, = 2F,,,): 


1+ 


for Y > 2, the error, with an accuracy of 10%, is determined from the formula 


9.2 \¥ 

pea = — + SS)”. (16) 
As is clear from (16) and (13), a shift of the signal frequency band to the low-frequency region decreases 

the error, For latge y , the error is decreased by a factor of approximately V2, In practice, however, it is not 

always possible to use just the low-frequency portion of the spectrum. 


Previously we considered transmission by triangular video pulses, It may be shown, by analogous methods, 
that with the same dynamic range, the error for PFMI-AM, i.e., for transmission by radio pulses with a triangular 


envelope of amplitude U,,, for y > 2, will be 
2 


B. Pulse- frequency modulation with constant porosity (PFMIL). 


With PFMII, the width of the pulses and the frequency with which they follow one another, vary propor- 
tionally to the magnitude of the measured parameter in such a manner that the porosity remains constant and 
equal to two, 


The signal is 


AQ, th= (t—te), 


= 
PY 
| 
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where F (t,) is a function describing the individual pulses and depending on the parameter A, 


The square of the derivative 9@AQ,t)/ @ can be determined by the same method used for PFMI, 
Computations show that. § pry = 5 


Table 1 gives the formulas for computing the relative mean square error for various cases of PFM, 


TABLE 1 


Form of Mean-square-error for 
modulation weak noise 


0.25 


AM 0.25 
Vy 


FM 
FM —AM 


PF MI 9.2 

fi = 2Fm< fe, Y>21) oy 
PF MI. 

hi. im 


PFMI-AM 


PF MII 
fi = 2Fm< fg, >2 


PF MI ‘ 10-°.8 -0.9 Bp 
hi fm >fe 3.5< Bp <25 


Note, The text provides the limits of applicability of these formulas, 


Based on the results obtained, the following conclusions may be drawn. 


1. For the same dynamic range, PFMI is equivalent in noise stability, to PFMIL, In noise stability, 
PFMI-AM is worse than PFMI by a factor of 4 ¥2, 


2. In all the cases of PFM considered, the magnitude of the error depends in the same ways on p and on 
Y, where this dependence is the same as with FM, 


3. For operation in the low-frequency portion of the spectrum (f,=2F mp « fy), the noise stability of PFM 
is somewhat higher than that of FM, 


. For operation in the high-frequency portion of the spectrum fj, fp, > fg, PFM is equivalent, from the 
point of view of noise stability, to FM. 


Pulse-Time and Pulse- Width Modulation (PTM and PWM), 


With PTM, the signal is a sequence of pulses of constant form which are shifted in time with respect to the 
interval markers, these shifts being proportional to the parameter A. Let the period of the marker points be T 
and for variations of the parameter’ from — 1 to +1, let the time shift of the pulses vary from — T/2 to 
T/2, Denoting the pulse envelope by U(t), we write the signal in the form 


A(a, t)= SU +57). 


k=0 


0.4 ‘ 
(y—1)’ PTM-AM — 
| 
e(y—1) for 0.35 
4 PWMI. 
log p> 
| 0.49 
PWMI 4.4 
1.7 3 9.2\!/ 
i 
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To determine the derivative 9U/0, we carry out a change of variables. t' =t— (kT + TA/2), Then, 


2A.) | 7 


For simplicity of analysis, we consider transmission by triangular pulses with amplitude 2U,, and width r. 
At the moment of pulse occurrence, 


(ur+ Za) 

4U2,T? 

(xr+ a)-= 


Substituting the magnitude of Ip; y in the expression for the error, we obtain 
PTM Y32U,,T 

For PTM, the parameter y will be 


0.5 
$rm (17) 


Expressing the error with PTM in generalized parameters, we get 


Computations show that transmission by radio pulses, i.e., by the use of a secondary AM on the carrier 
frequency, will entail an error, for the same dynamic range and band, four times as great: 


Thus, exactly as for FM, the error in PTM is determined by the parameters p and y. However, unlike 
FM, the error in PTM varies in inverse proportion to the square root of the frequency band- width taken by the 
signal, With PTM, the error is completely determined by, the steepness and width of the pulse edge. It remains 


unchanged if transmission is by pulses with a trapezoidal envelope having the same edge as the triangular pulses, 
since the derivative AQ, t)/@A equals zero on the flat portion of the pulses, 


It is easily shown that, with one-sided PWM the error will be larger by a factor of ¥2, then with PTM, if 
the pulses have the same steepness and width of edge: 


0,35. 


For two-sided PWM, the error will be larger. by a factor of ¥2, than for one-sided PWM: 


Correspondingly, 
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Pulse-Code Modulation (PCM) 


We now estimate the noise stability of parameter transmission by a simple binary code which uses all 
possible combinations, The errors in such transmission arise both from quantification and from the suppression 
of correct, or generation of incorrect, pulses. We consider errors due to misrepresentation of the code pulses 
by noise, Let the code consist of pulses, the values of which increase with increasing ordinal number of the 
pulses in accordance with the binary law: 2, 2', 2*,,.., 2%" *, 29-4, We assume that the probability that 
one pulse is misrepresented by noise is py, where suppression of a correct pulse and generation of a false pulse 
are equiprobable, Then, the probability that just one pulse is misrepresented in a code combination consisting 
of n pulses equals p(1)=np, (1— p,)"~4, For np, « 1, p(1)=np,. 


We assume that the probability of two misrepresentations occurring in one code combination satisfies the 
inequality p(2) « p(1), We now determine the mean square error under the condition that misrepresentation 
is equally probable for any pulse in the code. The mean value of the squared error, if one pulse in the code 
combination is misrepresented, equals 


By taking into account that the entire parameter range, using pulse-code modulation, is transmitted by 
means of (2"— 1) gradations, and that the probability of one error is p(1), we determine the mean square 
relative error: 


_ 


1 


(1) Pi =A(n) Py (19) 


For n>» 1, A(n) = 


The probability that one pulse of a binary code be misrepresented is determined, as shown in [1], by its 
energy, the method of transmission and by the specific noise power, For any noise level, 


where Q? is the maximum pulse energy equal, when the envelope is rectangular with pulse duration r, to 


e- 
and y is a coefficient which depends on the method of transmitting the binary digits (form of keying modulation) 
(3), 
The function V is defined by the expression 


V (z) = me 


1 
x 


dz. 


For amplitude modulation, v ,y4; = 1, for frequency modulation, vpyy, = ¥ 2 and for phase modulation, 
’pmt = 2. We will assume that the code pulses are transmitted with porosity 


T 


bed at’ 
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1.4 1.96 
4"—1 
where A (n) = — 


We introduce an additional parameter of code transmission, 


v 
and write 


Substituting (20) in (19), we get 
A(nyv (2), (22) 


It follows from (21) and (18) that, for a given value of the generalized transmission parameter p, which 
does not depend on the method of transmission, the error due to noise in code modulation increases with an in- 
creasing number of pulses in the code, We now determine the error due to the quantification of the signal, We 

4m 4m 


assume that all values of this error are equally probable in the range — ‘teen , where 


1 
bn = (23) 


Then, the probability density of the error equals —1_., and the mean square value of the error is deter- 


4m 
mined from the formula: 


By the use of (23), we obtain 


As was to be expected, the error due to quantification decreases as the number of pulses in the code in- 
creases, If we assume independence between the misrepresentations due to noise and due to quantification, we 
obtain the total mean square error by addition: 


1 


2 


Be 
4+ A(n)V (2) (24) 
Since the quantification error decreases with an increasing number of code pulses and the noise error in- 
creases with increasing n then, obviously, there is some optimal number of pulses in the code for which the 
total error is a minimum, An analytic solution of the equation, a5bcom / dn =0 is mathematically difficult. 
For the range of values, 3.5 < 6p < 25, the functional dependence, nop = f (6p), was determined by computa- 
tion, With an accuracy of 7%, this dependence, in the range given, can be approximated by a linear function: 


= 045 + 0.6 Bp. (25) 


(21) 
4m Am 
2 { 33 2 A? 
v= | =F. 
_4m 
2 2 
| 154 


With an accuracy of 5%, the minimum value of the error obtained for n=MNop: in the range of Bp given 
above, can be calculated from the formula 


= (26) 


On the basis of the results obtained, the conclusion may be drawn that the noise stability of PCM is com- 
pletely determined by the generalized signal parameter p, and code transmission parameter §. The greatest 
value of 8, and, consequently, the least error, is guaranteed by phase keying without intervals (the magnitude 
of 8 depends on the method of keying and on the code porosity), For each value of the product gp, there is an 
optimum number of pulses in the code, for which the error will be a minimum. It is obvious that, for a given 
length of code, the frequency band taken by the signal is proportional to the number of pulses in the code: 


1 una 
Afecm= = (27) 


where p is a coefficient characterizing the dependence between the duration of transmission for a binary symbol 
and the frequency band necessary for its transmission. For PCM the parameter y will be 


= (28) 
The optimal number of pulses in the code corresponds to the optimal value of the parameter y : 
YpcM(opt) = (0.9 + 0.6Bp). (29) 
Figure 2 shows the dependence of the minimum total error with PCM on the magnitude of p for the cases 


of amplitude keying with porosities of q, = 2 and &p =1, frequency modulation without intervals (ap = 1) and 
phase modulation without intervals (ap ='). 


Figure 3 gives the dependence of Yopt on the magnitude of p for various methods of keying. In com- 
puting Y, it was assumed that pant = 1-3 and perme = 2.6. 


Comparative Noise Stability for Various Transmission Methods 


The table gives the definitive expressions for the error in the presence of weak noise, It is clear from 
the table that FM, PFM, PTM and PWM permit a decrease in error, as compared to AM, at the cost of broaden- 
ing the frequency band used, FM and PFM are more 
effective than PTM and PWM since the former allow a 
40 60 oO 
gas Yt decrease in error inversely proportional to the first power 
Pair; of the frequency band- width whereas PTM and PWM allow 
a decrease in error which is inversely proportional to the 
square root of the frequency band-width, These results 
differ somewhat from the conclusions of work [1), where 
it is shown that PTM and FM are equivalent from the 
| point of view of noise stability when there is fluctuating 


\ noise, The discrepancy in results arises from different 
\ 6 conditions of comparison, 


In work [1], the modulation comparisons are made 
Fig. 2. Dependence of the mean square error for equal signal powers, while in this work, the condition 
on the parameter p for: 1) AM; 2) PCM-AMt, is that the dynamic ranges be equal, Systems with pulse- 
a, = 2; 3) PCM-AMt, a, =1; 4) PCM-FMt, frequency modulation, in many cases, permit a somewhat 
a, =1; 5) PCM-PMt, ap=1; 6) FM with higher degree of noise stability than systems with frequency 
optimal band; 7) FM for the band optimal modulation, although this increase is small, Therefore, in 
(for the same p) for PCM-AMt for a, = 2; 8) practice, it may be considered that, with weak noise, 
FM for the band optimal for PCM-AMt for systems with pulse-frequency modulation are equivalent 
Op =1. to systems with frequency modulation, from the point of 

view of noise stability, 


We now consider in more detail systems with frequency modulation, and compare them with systems with 
code modulation, In systems with code modulation, as was shown above, for each ratio of signal level, noise 
level and speed of action, there is a certain optimal value of number of pulses in the code and corresponding 

to this, an optimal value of signal band- width, for which mini- 
¥ opt mum error is attained, It follows from (8) that, for systems 
190 with frequency modulation, the error can be made arbitrarily 
small by broadening the frequency band (by increasing 7 ). 


2 However, this is not entirely the case. Formula (8) is valid 


for weak noise, When the band is broadened the noise in- 
creases and, for a definite value of y , Formula (8) becomes 
invalid, In work [4], correlation methods were used to investi- 
Fae gate the noise stability of an idealized frequency receiver, 
A consisting of an input filter, an ideal frequency discriminator 
Ay a and an output filter, The mean square error with parameter 
transmission and reception by this receiver were determined 


U; 
7) Wp > 5and > 2.2 
Fig. 3. Dependence of the optimal value » 
of the parameter y on p for: 1) FM; by the formula 
2) PCM-FMt(q@,=1); 3) PCM-PMt 
=1); 4) PCM-AMt =2); 5) 2 _ 0.083 0.26 ay 
AMt (cp = 1). (30) 
< 10-1+1-8 log (31) 
then the second term in (30) is much less than the first term, and 
0.3 


Expression (32) almost coincides with (8), and characterizes the error for weak noise, Inequality (31) defines the 
limit of applicability of the formulas for weak noise [Formulas (8) and (32)). For 3 < y < 200, this limit corre- 
sponds to a signal-to-noise ratio at the receiver input of (U;/Up);, > 3- 3.5, The second term of Formula (30) 
becomes the dominant one for weak noise, . It follows from (30) that for each value of p, there is an optimal 
value of the parameter 7 , and a corresponding frequency band, with which the minimum value of error is 
obtained: 


The minimum error, for y = opt» can be determined from Formula (30), 


Thus, it is possible to note certain analogies between the properties of code modulation and frequency 
modulation, In both cases there is an optimal signal band which guarantees a minimum error when transmitting 
over a channel with noise fluctuations, It is interesting to compare the minimum values of error for FM and 
PCM, Figure 2 shows the dependence of the minimum error on the parameter p for PCM for various secondary 
modulations and for different code porosities, as well as the dependence of the minimum error on p for FM, It 
is clear from the Figure that for 7 <p < 18, transmission by FM provides a lower error than transmission by 
PCM with the most stable, under noise, secondary keying, namely, phase keying (PCM-PMt), For 7 < p < 31, 
frequency modulation gives a lower error than pulse-code modulation followed by frequency keying (PCM-FMt?), 
and for 7 < p < 60, FM is better in this sense than PCM followed by amplitude keying (PCM-AMt), On the . 
basis of these results, the following general conclusion can be drawn, that for very weak noise (§ < 107° - 10°), 
pulse-code modulation has higher noise stability than frequency modulation, In the range of error, 10*< §< 
< 107° - 10°, frequency modulation permits a lower transmission error to be obtained than does pulse- code 
modulation, 
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We now compare FM and PCM from the point of view of band- width taken by the signal, Figure 3 shows 
the dependence of the parameter y opt = Afopt/2F, on p for FM and PCM, It is clear from the figure that 
in the range of values of p for which FM is better than PCM, frequency modulation requires a narrower band- 
width than does pulse-code modulation, By comparison with PCM-AMt, frequency modulation for the optimal 
band requires, for p > 22, a larger band-width, It is interesting to compare the noise stability of FM and PCM-AMt 
for equal band- widths. Figure 2 shows the results of error computations in systems with FM for the frequency band 
equal to the optimal band for PCM-AMt, with Op = 1 and a, = 2, The curves show that, for identical frequency 
bands, FM is more stable under noise than is PCM-AMt (ap =1) for 7 < p < 50, and PCM-AMt (ap = 2) for 
1< p< 80. 


Thus, for each form of simple binary code, there is a certain range of values of p in which systems with 
frequency modulation are more stable under noise and require narrower band- widths, 


We consider an example, Let it be given that U, = 50 millivolts, o = = and T = 0.1 second, 


It is required to determine the mean square error and the corresponding signal frequency bands for AM, FM and 
PCM, We find the magnitude of generalized parameter p from Formula (4), Since Fm, = 1/2T, p =11.2. 


According to (5), the error for AM will equal § 4), = 12.6%, 


The frequency band for AM equals Af,y, = 2Fm = 10 cycles, With this band- width, the signal-to-noise 
ratio at the receiver input will equal 


To compute the noise stability with FM, we determine the optimal value of the parameter y and the 


optimal deviation frequency, According to (33), y opt = 7.8. 


We take y =8 and, correspondingly, f, = 40 cycles, With this, the frequency band- width taken by the 
signal will be 


From (30) we determine the magnitude of the mean square error with FM: 


2. = 20.2 x 
20.2 x 10-* 


From whence, 
= 0.45%. 


With a 90-cycle band- width, the signal-to-noise ratio at the FM-receiver input will be 


(z*) 
U EM 4f 


We now determine the noise stability of PCM with secondary amplitude modulation (i,e., by the presence 
or absence of a code pulse) and with the intervals between pulses the same length as the pulses, With this, y = 
=1, = 2 and g, determined from Formula (21), equals 0,25. 


The product Bp , equals 2.8. Consequently, it is not possible to use Formulas (25) and (26) for determining 
the optimal number of pulses in the code, or the error, By computation using Formula (24), it is easily found 
that Noy, = 3, With this, the squared error equals 6 bcyy = 0.0017 + 0,428V (1.62). 


From tables,* we find that V(1,62) =0.053. Hence, 6 pcyy- amt = 15.6%. The band- width required for 


2? 
. \ e 2d: are given, for example, in the book of B, V, Gnedenko, “Course 


* Tables of the function  (z)= 7 
in Probability Theory," [In Russian]. 
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PCM-AMt for y =1.3 [see (27)] is Afpcm- AMt = 156 cycles. The signal-to-noise ratio at the PCM- AMt -re- 
ceiver will equal 


U 
3 = 1.6. 


4fpcm-AMt 


(U, PCM-AMt = 


We now determine the noise stability of PCM with secondary frequency keying and without intervals, i.e,, 
fory = ¥2 and 


In this case, 8 = 0.5,8p =5.6 and, for the determination of Dopt and §pcM-FMt » One can use Formulas 
(25) and (26): Dopt ws 4 and 5 PCM-FMt = 4%. Band- width, AfpcmM-FMt = 208 cycles, At the input of the 
PCM-FMt receiver, the signal-to-noise ratio is 


(U5 / Un kin POM-FMt = 


Thus, with the data considered, FM gives better noise stability than PCM-AMt and PCM-FMt, and requires 
& narrower band-width, In comparison with amplitude modulation, FM requires, for the optimal deviation, 9 
times the band- width, but gives an error which is less by a factor of 29. 


It sometimes happens in practice that, in FM systems with very low speeds of action (on the order of a 
second), it becomes difficult to attain the optimal band- width due to frequency instability of the transmitter 
and the frequency characteristics of the receiver, and deviation frequencies greater than optimal occur. In such 
& case, an FM system may turn out to have less noise stability than a PCM system, in which it is generally 
easier to attain an optimal band-width. However, in such slowly acting systems, the error due to noise is generally 
very small, with the signal-to-noise ratio at the receiver input greater than two, so that the requirements can also 
be satisfied by a system with y > y,,;. Thus, for example, in an FM system with fy = 75 cycles, T = 0.5 seconds 
and (U,/U,);,, = 2.8, the error is 0.24%, With this, y = 75, The optimal value of y for the given system is 40, 
1.€., fy op¢ #40 cycles, Computations show that the error for Y = __, and the same values of signal, specific 
noise ond epend of action, would equal 0.04%, i.e., would be six times smaller, If PCM-FMt (vy = ¥2, a, =1) 
were used for parameter transmission, then the mean square error for the same data would amount to 0.06%, 
i,¢., would be somewhat larger than for FM with f = 40 cycles, but would be four times less than for FM with 
fq = 75 cycles, 


In this work, noise stability for PFM, PTM and PWM was considered only for weak noise. This hindered 
comparisons of them with other types of modulation, The results of analysis of noise stability of PTM, PWM 
and PFM with relatively strong noise fluctuations will be published in subsequent papers, In these will be given 
the limits of applicability of the formulas for weak noise for these forms of modulation, and comparisons will 
be made of PTM, PWM and PFM with FM and PCM, 
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LINEAR TRANSFORMATIONS OF BINARY CODES 


N, Ila, Matiukhin 


(Moscow) 


Linear transformations of binary codes with fixed number of elements per code 
are considered such, that there is a one-to-one correspondence between any symbol 
and an object of coding (nondegenerate codes).* 


Parallel to the binary number system with its wide field of applications in digital electronic computers, 
telemetry and telecontrol, a large number of other binary codes with their own distinct characteristics is also 
known, Here belong Hamming's self-correcting code [1], Gray code [2] for converters of continuous quantities 
into discontinuous ones, or the channel capacity increasing Fano-Shannon code [3], etc, 


The methods of separately obtaining each of these single codes,as well as their properties,are described 
in detail in literature, However, much less attention has been paid to methods of construction applicable to a 
wider class of codes, 


One of these methods based on representing coding symbols as vectors in a multi-dimensional space, and 
on investigations of linear transformations in that space is considered below, 


We shall assume that numbers are the only objects of coding. To every number M from the set of numbers 
permissible in a given system there corresponds an encoding symbol of m elements 4, Og, ... » Sp» (m-tuple), 
where each element ofthe symbol a; can take only one of two values (0 or 1), The total number of different numbers 
which can be obtained using a binary m-tuple is equal to the number of different encoding symbols: 


where Ci is the number of selections of i elements out of m, 


The number of different correspondences (codes) between 2™ numbers and 2™ m-tuple code symbols is 
equal to the number of all possible rearrangements of symbols N,, = 2™!.,When m = 4, they already number 10", 
approximately, 


We are choosing the 2™ numbers M to be the first integers, starting from zero: 0,1, 2,..., 27-1. 


To any number M shall correspond m-tuple (04, a, . . « » &p) by M = >} «2°, that is by making 


i=1 
use of the binary representation of number M, 
The code thus obtained is called normal code Xp. 


* Self-correcting codes, for example, do not belong to this group, as in this case one object of coding corresponds 
to several encoding symbols, 


Any other possible code X can be written in the form of a substitution table establishing a correspondence 
between the normal and the other code, for example when m = 2 


01, 10, 11 On 01, 10, 44 , ingeneral Xn 
10, 11, 00, 01 11, 10, 01, 00 X/}* 


Linear Substitutions Subgroup and Its Characteristics 


We shall study some properties of the interesting substitution group. To do this we shall consider the 
m-dimensional vector space over field P;, the residues modulo 2[4), In this space any vector is in the form 
X = (xq, Xg,  » + » Xp), Where x; is equal either to 0 or 1, The addition of two vectors and the multiplication 
of a vector by a number from field P, is given by formulas 


AX = (Ady, » Am). 


Null vector 0, such that x +0 =x, is of the form 0 =(0, 0, ... , 0); vector x" opposite tox, i.e., such 
thatx+x‘'=0 is equal to vectorx itself (x +x=0). 


The following m linearly independent vectors are selected in our space: 
e, = (1, 0, 0, eee , 0, 0), 
e, = (0, 1, 0), 
e, = (0, 0, 4,...,0, 0), 


and the set of these vectors is called the base of our vector space. Any vector x can be obtained as a linear 
combination of the base vectors: 


In the m-dimensional space there are 2™ different vectors and therefore a code symbol of a binary code 
may be made to correspond to each vector, We agree to write the symbol corresponding to a given vector in 
such a way that the i-th element of the symbol from the left corresponds at the same time to the i-th component 
of the vector and to the i-th coordinate in the base €;, €3,...,€,). 


For briefness we shall also denote a vector by its code symbol remembering it, however, when performing 
Operations on vectors, For example, instead x =(1,0,0,1) we shall write simply x = 1001, 


The vectors of our space can be subjected to linear-transformation represented by square matrices of order 
m_, having all elements either zero or unity, 


The operation of linear transformation is denoted in the form y = Ax, where x is the input vector, andy 
the vector corresponding to it via transformation given by matrix A. 


Transformations are performed in accordance with the usual matrix by vector multiplication rule, For 
example, if 


1114 
A=|101|, 
100 


then y=Ax=1-11-0@1-1, 1-1@0-0@1-4, 1-14@0-060-4=001. 


* The sign @ means addition mod 2, 


em = (0, 0, 0... 40, 4), 
| 


The determinant of matrix A is found by applying usual rules with one modification viz, in any formula 
for the evaluation of the determinant all negative signs are changed to positive,* 


Multiplication and addition of matrices are also performec according to known rules, and have properties 
analogous to those of common matrices, 


If the determinant of matrix A is not zero (matrix A nondegenerate), there exists an inverse matrix A™* 
such that the product A- A“! =E, where E is the unitary matrix, Later, we shall be mainly interested in non- 
degenerate matrices, 


We shall consider the following two examples of nondegenerate transformations, 


1. Identical transformation Ex =x 


1000 
0100 
E=l9010 
0001 


2, The interchange of rows transformation 


1000 
0001 
0010 
0100 


which interchanges the second with the fourth row, It is obvious that A“ A, i.e., A’ =E, 


How the nondegenerate transformations determine the one-to-one correspondence between the input 
vector x and the transformed vectory can be seen by regarding the vectorsx andy as code symbols and 
using A to form a substitution table: 


000, O01, 010, Off,... 
Pes A(001), A(010), 


AX. a? where A is any nondegenerate matrix, are fully determined by their 
matrix A, Such substitutions shall be called linear substitutions, 


Any code obtainable by linear substitution shall be-called linearly dependent on the normal code, The 
number of such codes (including the normal code) is equal to the number of different matrices A: 


Substitutions of the type ( 


m(m—1) ™ 


Na =2 (2'—1). 
i—1 


The above formula is derived in the Appendix. The number of different matrices when m = 4 is already 
-104, 


It has been shown above that all the possible 2™ substitutions form a finite group, One can show that all 
the linear substitutions form a subgroup of that group as all the possible nondegenerate matrices form a group 
closed with respect to matrix multiplication, The substitution group can therefore be decomposed into N= 
= 2my/ Nv? (right-hand or left-hand) cosets relative to that sub-group [5], A coset is determined by one of its 
codes int the substitution table), All other codes of that coset may be obtained by linear substitutions carried 
out on that code, 


* Subtraction in P, field is equivalent to addition as x = ~ x, 
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Let us have, for example, a substitution ( x 


) , hot belonging to the linear substitution sub-group. By 


xX. 
means of this substitution a coset ~ ) can be constructed, The linear sub-group alsoforms a coset generated 


by identical substitution ae Therefore, any nondegenerate code may be formed by one of N codes (sub- 


stitutions tables) defining N different cosets, and by one of N‘;' matrices from the linear sub-group. In this work, 
however, as mentioned above, our consideration is limited only to codes of the linear sub-group, 


Later, it will be necessary to form matrix A from a linear substitution table, Matrix A is fully given by 
the simultaneous equations: 
Ye = Ax,, 


Yn = AXm; 


where vectors and also yy, entering the substitution table in twos, are linearly 
independent, Otherwise either matrix A would be degenerate or not determined,due to an insufficient number 
of equations, 


If instead of vectorsX4, Xz, ...+%Xm we take normal base vectors Cg, . » then vectors yy, 
identical with columns of matrix A, If 


Examples of Encoding Systems of the LinearSub-group and Physical 


Realization of Linear Transformation Matrices 


Example 1, We shall construct transformation Aq such that the even symbol of code X, and code Y be 
identical and that the odd X, symbol have in Y a symbol obtained as a complement of its predecessor. The 
X, symbols are listed according to Formula (1), When m = 4, our rules can be formulated in the following way: 


ga = Yont1 = Xan 1111 


A, 43) = 1111. 


But Xy_ @ Xyn4, = 0001 (they differ only in the least significant digit), that is 


A, (0004) = 1111. 


then 
(2) 
or 
162 


One arranges the other three conditions for Ay as follows: 


ail A, (0010) == 0010, A, (0100) = 0100, A, (1000) = 1000. 
It is not difficult to verify that these four rules satisfy the requirement for vector linear independence, as 
stated above, 
2 Using Relations (2), we obtain matrix 
rk, 1001 
0101 
oot 
0001 
and the corresponding table 
0000, 0001, 0010, 0011, 0100, O101,..., 4190, 1111 
0000, 1111, 1101, Of100, 1011,..., 1110, 
Example 2, We shall construct transformation A, such that to the adjacent symbols x,, and X44 of code 
X,, shall correspond y, and yp,, of code Y, differing from one another in one place only, First we shall find 
how symbols x,, and x,,,; differ from one another, in other words, what is their sum Xp ® Xn44. 


tly 
| Symbol X41 is an arithmetical sum of x, and the symbol which has 1 only in the least significant place, 


When adding, the carry-over 1 may have to be carried several places to the left and, depending on how many 
places it has been carried, the sum of interesting symbols may assume one of the following values: 


. 0001, 
Xn 
1111. 
As assumed above, 
Yn ® 


where e; is the symbol which has unity only in the i-th place, that is a normal base vector. 
In view of As (x, ® Xn41) =Yn Yay, We can write down the system of equations for A; (when m = 4): 


Az (0001) = 0001, A, (0011) =v010, A_(O111)= 0400, Ay (1111) — 1000. (3) 
From these equations the matrix of the inverse transformation Aj can be determined straightaway, 
writing these as follows: 
0001 = Az1(0001), 0011 = (0010), O111 = (0100), 4441 = AZ*(4000). 
Hence, in agreement with (2), we have 


1000 


_, _ |1100 
Ay = | 4410}: 


y? 1144 


To find Ag it is convenient to form another system of equations by adding the equations in (3) in two's 
in such a way that the base vector transforms appear on the left-hand side, In this way, we obtain 


Az (0001) = 0001, Ay (0010) = 0011, A, (0100) = 0110, (1000) = 1100. 


Hence, 
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4 
1000 
1100 
=| 
= 


The corresponding substitution table takes the form 


0000, 0001, 0010, 0011, 0100, O101, 0110, O111,..., 1441 
0001, O011, 0010, 0110, 041141, 01014, 0100,..., 1000) 


The code thus obtained is called Gray code and is used with a quantizer transforming continuous quantities 


into discontinuous ones by using a code wheel or a plate, 


Fig, 1 


Figure 1 shows a schematic outline of a shift-digit 
converter with a Gray code wheel, The angular position 
of the shaft is converted into a number with an accuracy 
of 1/16th of a revolution, The dotted places on the 
whee] are electric conductors, The number correspond- 
ing to the position of the wheel is brought down from 
contact 1-4 in the form of a binary digit. Voltage is 
brought to the wheel through a central contact, On 
each contact 1-4 voltage corresponds to digit 1, no 
voltage to digit 0. 


While the voltage in the contacts is changing 
the error in the determination of the wheel position, 
which may occur in the normal code, is eliminated 
because the adjacent positions on the wheel differ in 


one place only, For instance, normal code symbols corresponding to numbers 7 and 8 differ in all their places 
and during transition from position 7 to position 8 any symbol corresponding to one of the positions 1-16 might 
be obtained, depending on the diameters of the contacts 1-4 and how well they were fixed, 


Figure 2 shows a photograph of a wheel with 512 positions, obtained photographically, 


Fig, 2 


It should be noted that though it follows from 
(3) that there exist m! different transformations Ag, 
which can be obtained by rearranging the right-hand 
side of (3) in different ways, the codes thus obtained 
differ from Gray code only by transposition of places 
and therefore are of little interest to us, 


The inverse transformation, from Gray code into 
the normal code performed by matrix A3*, can have 
a very simple technical realization, 


The converter (Fig, 3) has trigger T fed through 
its counter successively — place after place — by Gray 
code symbols in the form of series of pulses, starting 
from the top digit, The sum of digits (modulo 2) of 
the higher places of Gray code symbol appears con- 
secutively on the trigger. The sum of the k highest 
places of Gray code symbol gives the k-th place of 
the normal code symbol if our adding starts with the 


top digit, The normal code symbol is being obtained place after place from trigger T through valve K, It is 
easy to verify that this code transformation corresponds to matrix A3’, 


The demerit of the converter is that the symbols are obtained starting from the top digits on the left of 
the binary representation, and that these symbols cannot therefore be used directly in automatic digital com- 
puters where the feeding of the code symbol always starts with the least significant digit on the right of the 


binary representation, 


In the example below properties will be considered of a code whose normal code symbols are obtained 


by using a device shown in Fig. 3, which begins with the 


Example 3, We shall examine the properties of a 
form 
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least significant digit of the binary system, 
code whose inverse transform matrix is given in the 


~2 
K-04 =~ 
Sue 
aw’ 4 
10 
7 
|| 
| 


on 
cy 


1111 


normal | Gray code Y 0 
code X 
_ cyclic But matrix Aj? is the transpose of Aj! matrix 


in our last example, hence matrix Ag is obtained as 


a transpose of Ag matrix of the same example: 


pulses 
A AA A 
pulses | | | 0110 
| A= 
coy [AL i Ai 
7 
rT O | | | Pi Transformation Ag leads to relations: 
A A As (0001) = 0014, As (0011) = 0401, 
code 1 Ay (0114) = 1001, Ay (1144) — 0001. (4) 
Fig. 3 


In this way any two consecutive y,, and yn44 
are distinguished by the least significant digit and 
differ in not more than one place, Comparing (4) 
with (3) we find that the symbols of system Y have 


8 A; — with the exception of the least significant place ~ 
code the same structure as Gray code symbols, 
Fig. 4 Example 4, We shall describe a method of 


converting Gray code symbols into normal by using, 
as an intermediate device,the converter shown in 
Fig. 3 which is fed with Gray code symbols starting 
" with the least significant digit, The converter in 
this case performs the transformation corresponding to the matrix in Example 3: 


1444 

_ | 0111 
=| 
oovt 


We shall find an additional transformation B such that when applied after Ag* (Fig. 4) the normal code 
will result, 


As =A3'y and x =Bsg, then x =BAj! y 

Hence BAs!=A3! and that is when m=4 
1000) 11100 1100 
1100| |0110 1010 


(5) 
114114] |o004 1000 


It can be seen from (5) that transformation B is such that top digit z, being at the same time the least 
significant digit xg, forms the digit of the required symbols 2, Zg, Zz when added successively to z, modulo 2, 
The realization of this transformation is shown in Fig. 5. 


We shall now consider relay contact networks suitable for obtaining a realization of any transformation 


matrix (Fig. 6), In the case of active zero action each matrix row represents one switch circuit which gets voltage 


- 

| 
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on its left input terminals representing digits 1 and 0, respectively. On the right terminals the representation 
of the digit yy, of the transformed symbol y is obtained; if on the terminal y, voltage represents 1 and on the 
terminal yy, it represents 0 then Yk =1. The reciprocal combination of voltages means that y, =0. Should the 


code Z 


» 


normal AA 


Fig. 5, DID — disaccord decoder, DI, — delay line for (n+ 1) 
places. 


44 


+ 


> @ Y, 0 xX xX 


% 
4 % % 
100 
Get Fig. 7. = 110 z= 
4, 
7,20 = 
0 digit be represented by the breaking of the circuit, the 
4014 lower left and the upper right terminals can be left out. 
Fig. 6 Aq | 0110 The i-th row switch circuit is expressed by successive 
nageee nosy couplings of the contact circuit corresponding to the uni- 


tary elements in this matrix row. The contacts of any 
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such circuit are controlled by the k-th digit of the code being transformed, The network diagram and the usual 
terms are given in Fig, 5, 


It is ‘not difficult to ascertain that successive couplings of such circuits form a network which adds modulo 2 
tothe digits xy, X2,..., Xk, by which the contacts are controlled, 


The general scheme shown in Fig. 6 can be simplified considerably when the columns of matrix A contain 
a large number of 1's, It can be best seen how it is done in the example of the construction of the matrix which 
is inverse to Gray code transformation, In Fig, 7 the general and simplified schemes of this converter are given, 
The least number of elements required for any nondegenerate matrix cannot be less than its order m, 


SUPPLEMENTARY REMARKS 


One of the important aims of the coding theory is the construction of codes on which various operations 
will be most easily performed (for example, arithmetic operations, etc.) as well as the construction of algerithms 
which will perform these operations, 


It can be shown (see Appendix, Section I) that the operation which in the nortral code corresponds to an 
arithmetical addition cannot be made simpler in any other code of the linear sub-group of binary codes, In fact, 
as proved in Appendix, Section I, in any other code (of the linear sub-group) it is necessary to perform three 
extra linear transformations in addition to all the elementary operations necessary for adding in the normal code, 


As far as the self-correcting codes are concerned,an example can be cited of linear transformations which 
lead from the normal to the self-correcting code, For instance, the transformation 


1000 
0100 
0010 
1110 
0001 
1101 


converts a four-digit normal code symbol into a seven-digit Hamming code symbol which corrects a single error, 
The inverse problem, however, the determination of the normal code symbol cannot be solved by linear trans- 
formation methods as the above matrix has no reciprocal, The solution of this problem requires more involved 
methods, 


At the same time it is possible, for example, to construct a single error-detecting code, Let x be the 
symbol of such a code and y = Ax, We shall only use those x symbols which have x,=0, Requiring an error 
in one single place to produce one in x, place of the x symbol, we obtain (when m = 4) 


(0001) = 1001, A*(0010) = 1010, 


A~*(0400) = 1100, (1000) = 1000 


14414 1444 
0100 0100 
0001 0001 


Looking at the matrix it is easy to see that the code obtained is the one which is known to check a single 
error by the sum of digits. In a general case, it is possible to construct a code with required characteristics by 
the method of code linear transformations, only when one is able to present these characteristics in the form of 
a one-to-one correspondence between linearly independent vectors of a given code system on the one hand and 
the required code system on the other, 


It is obvious from the above that the number of such coding systems is large and increases quickly with m, 
APPENDIX 


1, We shall determine the number of different nondegenerate transformation matrices, by applying 
mathematical induction, 


Let ae different transformations be in the (m— 1) dimensional space, We shall compute the number 
of matrices Ary of order m of the following type: 


A, = 


Here A,,.4 is a nondegenerate matrix, order (m— 1), 


For any matrices Ay,- 4 and B the matrix Aj, is nondegenerate, because the columns of Am- 4 are 
linearly independent and so are, therefore, the first m— 1 columns of matrix A. Also the m-th column must be 
linearly independent of the others because there is no linear combination of columns of Am-, which could give 
zero elements in the m-th column (this would be against the condition of linear independence), It follows also 


that matrices A4,, formed by all matrices Ay). and all B matrices exhaust all possible nondegenerate transform 
matrices A,, which have the above m-th column. 


- It fs clear that the number of such matrices is 
where Ng is the number of different B matrices, But B can be represented in the binary system by a number with 
(m— 1) digits and hence Np = 


Thus, the number of A,*, matrices with specified m-th column has been found, As we are free to fill in 
the m-th column in any way, with just the one exception of all zeros, and as with each such column goes the 


same number of A$, matrices as obtained originally, the total number of different transformations Arp is there- 
fore 


= (2 —1) N = (2™ — 1) 
As Ni, = 1 when m =1, we obtain 


m(m—1) m 


™m 
* 
Il, Performing arithmetical addition with transformed symbols. 


Let the arithmetic addition operation x + Xz =X, be performed on normal code symbols x and x, and let 
these be such that xg has the same dimensionality. 


It is required to find such an operation on the mapped symbols y that y,0y2=ys3, where yy = Ax, y2 = AX, 
Ys AXy 


We know that 
Xi +%2=% (6) 


where u is the binary carry-over code symbol, Hence 


yo @ Au, (7) 


which means that to be able to describe the arithmetical addition with transformed codes it is necessary to de- 
termine the transformed carry-over symbol Au, We shall show that Au cannot be expressed by yy and y, unless 
the inverse mapping A“? is used, 


To do this we shall represent u asu =u, @ Up, where ug is the symbol of unconditional carry-over and 
u, the symbol of conditional carry-over, 


To compute ug and uj, it is necessary to introduce, apart from addition modulo 2, yet another operation 
on the symbols; this operation is the place-by-place or logical multiplication x =x,x3. This operation is not 
linear with respect to transformation A, that is 


A (X1X2) Ax, Axe. (8) 


The unconditionsl carry-over symbol is obviously ug = M (x4xq), where 


010. 00 
001 . 00 
000 . 10 
000. 01 
000 . 00 


_ is the shift mapping transferring the obtained carry-over into the next place, 


The conditional carry-over vector is 


where u, = M [u,_, (x: xs)jandu, = M [u, (x: 


The first term of the sum w, corresponds to the conditional carry-over resulting from the unconditional 
one occurring in the preceding place. The second and following terms represent conditional carry-overs, each 
obtained from its predecessor, as seen from the last two formulas, 


Thus 
u = M [(x:x2) M [(xix2) (21 22)] @ M {[M (x1 (x1 @ - -)- (9) 


Hence, it follows from Construction (9) and Inequality (8) that to express the transformed symbol Au (7) 
in terms of yy and yg, it is necessary to calculate first vector u itself which can be done when y, and y, are 
known by substituting x, = A™4yy and x, = A““y, into Formula (9), As seen above, if the result of arithmetical 
addition of code symbols yy and yg is to be obtained in the Y code ~ following remarks on computing vector 
u and Au and using Formulas (6) and (7) — the same intermediate operations must be performed: computation 
of vector u , two inverse transformations A~’, one A transformation and two additions modulo 2, 


Further, if the result of the addition — expressed in the normal code— can be used then it is not necessary 
to perform transformation A when using Formula (6), Comparing arithmetical addition operation in the normal 
code (6) and any other code (7) it can be seen that in the last case it is necessary to perform three additional 
linear transformations: and Au, 
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STATIC CHARACTERISTICS OF TOROIDAL CORES AS A FUNCTION 
OF THEIR GEOMETRIC DIMENSIONS 


M, A, Rozenblat 


(Moscow) 


The paper treats of the influence of the ratio of the inside and outside diameters of 
toroidal cores on their magnetic characteristics and, in particular, on the fundamental 
magnetization curves and the forms of the hysteresis loops, Formulas are given for com- 
puting these characteristics, It is shown that core geometry is most essential in cases 
where magnetic materials are employed which have rectangular hysteresis act or sharp 
corners on the fundamental magnetization curve, 


INTRODUCTION 


Toroidal cores have found extensive employment in magnetic amplifiers, saturated chokes, dc measuring 
transformers, magnetic relays, magnetic memory devices, and in many devices where it is necessary to make 
maximum use of the magnetic properties inherent in the material of which the cores are made, The advantages 
of toroidal cores are that they allow practical exclusion of the deleterious influence of air gaps and leakage, 
and permit a more uniform magnetization to be obtained than do other types of cores, However, even toroidal 
cores do not allow complete uniformity of magnetization over the entire core to be obtained, This is a result 
of the differing magnetic path lengths along the inside and outside perimeters, thanks to which the intensity of 
the field generated by the current in the windings varies across the core’s cross section, 


The lack of uniformity of magnetization along the toroidal core’s cross section frequently leads to a 
significant deterioration in its magnetic properties, compared to the properties of the material forming the core, 


The present paper is devoted to a consideration of the influence of the ratio of outside and inside dia~ 
meters of toroidal cores on their magnetization curves and the form of their hysteresis loops, We remark that 
this question was partially treated in work [1], The authors of the work cited limited themselves to explaining 
the influence of the core dimensions on the steepness of the descending arm of the limiting hysteresis loop in j 
cases where alloys with rectangular hysteresis loops were used, With this, only the central (linear) parts of the 
descending arm of the hysteresis loop were considered, Here, this question is considered more fully. 


Magnetic Characteristics of Cores 


In what follows, we shall denote by b andh the magnetic induction and field strength at a given point 
of a core, and by B and H, the average values of magnetic induction and field strength of the core, Corresponding 
with this, the relationship, b = f(h), is the magnetic characteristic (magnetization curve or hysteresis loop) of 
the material, and the functional relationship, B = F (H), is the magnetic characteristic of the core, 


The fundamental assumption made in this work is that the core is taken from a homogeneous body, for any 
point of which the relationship, b= f(h), is valid. 


We denote by Ro and Ry, respectively, the outer and inner radii of the core, and by a their ratio, a= | 
=Ro/R;. Then, the average induction of the core is 


R 
i 
where b is the magnitude of the induction at the point situated at a distance R from the geometric center of the 
core, 


Analogously, we find the average field strength of the core: 


If the current I, in the winding w, placed around the core, generates a magnetizing force F =0.4 1 wl, 


then 
(2) 
By substituting this value of h in the equation for H, we find that 
Fink,/R Fine 
Thus, the average core field strength equals the field strength at points lying on a circle with radius 
h~inkO/R} 
called the mean harmonic radius, 


__ We mention that, in practice, the core field strength is frequently computed by dividing the magnetizing 
force by the average core length: 


Taw ot (5) 
which gives a lower value for the average field strength than does Formula (3), 
We provide several values of the ratio H/H', computed from the formula 


H a+i 


a i.4 1.7 2 3 4 5 
#/H’ 1.01 1.02 1,04 1.41 1,16 1.24 


For values of a = 2, it is possible for engineering computations to calculate the average core field 
strength by the simpler Formula (5) without any essential error.* 


In the future in calculating H we will be using Formula (3), . 


* In work (1), the authors used Formula (5) even for higher values of « which led to implausible results, For 
example, it turned out that an increase of « led in some cases to an improvement in the core’s magnetic 
characteristics, which does not, in fact, occur. 
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Core Magnetization Curves 


The core magnetization curve can be computed if one has an analytic expression for the material magneti- 
zation curve b= f(h), or can be found graphically if one has the function b= f(h) in the form of an experi- 
mental curve, 


We first consider the case when the material magneti- 


oe 7 geese 7 zation curve is represented by the broken-line curve 0-1-2 
oe paeey (Fig. 1). We denote by B, the saturation induction and by 
ag VA H, the field strength for which the core saturates, giving us 
a6 
h 
as b= (6) 
B, fo 
a2 If there are no saturated portions in the core, i.e., if 
h < Hy, then by substituting the value of b from (6) in 
a Tj ; a 3 254M, Equation (1), we find the average core induction: 
Fig. 1 B 4 B 
H, o H, 
I 
For a: =H,, which, according to Formula (3), corresponds to an average value of the core field 
strength of 
H=H,=4H, =, 


the inside layer of the core, having radius Ro, is saturated, Therefore, for H > Hy, it is necessary to compute 
the core induction by the formula 


Ro 
4 B, 
B = B, +\ 
P 
where Rp is the radius of the circle dividing the core into saturated and unsaturated sections: 


= 


If we substitute Expression (2) for h, and take (3) into account, we find, after integrating, the following 
expression for B when H = Hy: 


B 


Thus, for H < Hy, the core magnetization curve coincides with the material magnetization curve, and B 
varies in proportion to H, For H > Hy, the core magnetization curve becomes nonlinear, and may differ signifi- 
cantly from the material magnetization curve, 


We present the values of Hy, computed from Formula (7), for which the linearity of the core magnetization 
curve is disrupted, 


Ax/Hs i 0.84 0.78 | 0,73 | 0.69 0.55 | 0.46 
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Figure 1 gives, in relative units, the core magnetization curves computed by Formula (8) for a =1.4, 2, 
and 4, If, for a =1.4, the influence of the inhomogeneity of core magnetization is small, for a = 2 and a=4 
this influence becomes significant, and engenders a noticeable shortening of the linear portion of the magneti- 
zation curve and a smoothing out of the sharp corners of the curve, We note that the relative influence of the 
magnitude of a on the magnetization curve does not depend on the magnitude of the magnetic permeability 
(i.e., on the ratio Bs /H,), 


The character of the influence of « on the magnetization curve allows the assumption to be made that 
for materials without sharp corners, the influence of « will be less than for the case considered in Fig. 1. 


It is frequently convenient to approximate a magnetization curve without sharp corners by an arc tangent 
curve, sometimes with a linear term: 


b=A arc tg Ch+Dh, * (9) 
By substituting this expression in (1), and taking (2) into account, we obtain 


B= \ dR hdR, 


The second term in the right member equals DH, To simplify the first integral, we introduce a new 
variable: 


2=CF/2nR. 
Then, the integral takes the following form (3): 


dR = — de = — [ine —F init 
Replacing F by H in accordance with Formula (3), we find, after setting in the limits of integration: 


t+ (cu 
Ving 1+(cH 


+ are tg CHS — — are tg CHS). 


gauss 
5000 
4000 
0 ai a? as 4% Hoersted 
Fig. 2 


As an example, we consider the magnetization curve b = f(h), shown on Fig, 2 and computed. by Formula 
(9) with the values A = 3190 gauss, C = 45/oersted and D = 2400 gauss/oersted. This curve, for values of h < 0.5 
oersted, corresponds with high accuracy to the magnetization curve of alloy 80NKhS with maximum permeability 


= ‘tan’, ‘sh’ = ‘sinh’, — Publisher. 
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H max = 138,000 gauss/oersted and coercive force H, = 0.019 oersted (see [2], p. 541 ). The dotted curve is for 
the core magnetization curve B=F(H), computed by Formula (10) with a =4, Even for such a large value of 
a there is-a relatively small worsening of the magnetization curve. 


In certain cases, the material magnetization curves are approximated by hyperbolic sines: 


1 b 
| h= «oF arc sh Ch. (11) 
Substituting this value of h in (1), replacing R by the new x =1/R, and integrating, we find 


an expression for the core magnetization curve in the form: 


‘CH (a —1) 
B = ACH + 
a CH (a — 1) 
ailing 
A CH (a —1) CH (a—1) 
+ are sh — are sh SE (12) 


Computation of the core and material magnetization curves B =F (H) and b = f(h), respectively, by 
Formulas (11) and (12) also showed a relatively small difference between them for a = 4, 


gauss 


2 6 6 8 W @ 2M, cerned 
Fig. 3 


If the relationship b= f(h) is in the form of an experimental curve, the core magnetization curve can 
be found by graphically integrating Equation (1). For this, the values of H should be given, and the magnitude 
of the field strength h should be calculated at several points of the core situated at various distances R from 
its center. For these h, the corresponding values of b are found from the curve b = f(h), and the relationship 
b = g{R) is constructed for the given H, The area included between the axis of abscissas, the vertical lines R = 
= Ry and R = Ro, and the curve b = @(R) equals (when the scales are taken into account) the integral in the right 


member of (1), By repeating such a construction for various values of H, we find the desired functional relation- 
ship B =F(H), 


As an example, Fig. 3 shows, by the solid line, the measured fundamental magnetization curve b= f(h), 
for electrotechnic steel E42, and Fig. 4 shows the curves b = {R/R,), computed by the method described above 
for various values of H for a toroid with a=5, After graphical integration, the core magnetization curve B= 
= F (H), shown by the dotted line on Fig. 3,was obtained, 


Formulas (9) and (11) do not always assure a sufficiently accurate approximation of the fundamental 
magnetization curve in the regions of weak fields, For many magneto-soft materials, the magnetic permeability 


in these regions increases in direct proportion to the magnitude of the magnetic induction, and may be represented 
as 


+ kb, (13) 
where jg is the initial permeability. 


B, gauss 
14000 
12000 i7 


Fig. 5 


As an illustration, Fig, 5 shows the experimental 
dependence yp = f(b) for alloy 79NM (curve 1}, 80NKhS 
(curve 2) and the straight lines corresponding to Formula 
(13), For 79NM we have 9 = 20,500 gauss /oersted and 
k = 34 /oersted, and for 80NKhS, = 30,000 gauss/ oersted 
and k = 45 /oersted, 


We now explain the influence of the quantity a = 
=Ro/R;, on the initial segment of the core magnetization 
curve for material for which Formula (13) is valid, 


Since 4 = — , Formula (13) can be put in the form: 


R F kF 
By expressing F in terms of the field strength H in accordance with (3), we find that 
alna—kH (a — 1) 
B= In Ina—kH (a—1) a4) 


_ The dotted line on Fig. 5 shows the initial portion of the curve » =B/H =F (B), computed from Formula 
(14) for a core made of alloy 80NKhS (15 = 30,000 gauss /oersted and k = 45/ oersted) with a=4, It is clear 
that for a> 1 there is an increased core permeability in weak fields, as compared with the permeability 
of the material, 


If we denote by hy, the maximum value of field strength for which Formula (13) remains valid then, ob- 
viously, we can use Formula (14) for values 


Fig. 4 


Core Hysteresis Loops 


The hysteresis loop of a toroidal core, account being taken of the nonuniformity of magnetization along 
the cross section, is also computed by means of Formula (1), 


We shall first consider materials with a so-called “rectangular” hysteresis loop, represented in somewhat 
idealized form in Fig. 6. The residual induction B, is taken to equal the saturation induction Bs, The ascend- 
ing arm of the hysteresis loop is described by the equation: 


b=p(h—H.), (15) 


b where 


a le We shall assume that the core's initial magnetic state 
corresponds to point N of Fig. 6. Therefore, forh < H,~ B,/p 
we have b= — B, and, forh = H, +B,/y , the induction is b =B,, 
The core begins to reverse its magnetic polarity, from state — B, 
to state B,, when the core’s average field strength H reaches the 
value: 


| B 
H, Ina. (16) 


Fig. 6 With this, in accordance with (1) and (15), the mean value of 
core induction equals: 


I 


where R, is the radius of the circle partitioning the core into saturated (—B,) and unsaturated parts: 


For H = Hg, where 
+2 
H, = Ing, (18) 


the portion of the core adjoining the inner perimeter (R = Ry) has completely reversed its magnetic polarity, and 
has an induction of b = Bs. 


For H = Hg, where 
Ro B,\alna 
(19) 


all portions of the core have begun to reverse polarity, i.e., for H > Hs, at any point of the core b> — Bs. 
Therefore, Formula (17) is only valid for values of H < H, and H < Hy, The further course of the functional 


F 

an H,—-* 

After integration, we get the following expression for B, for H = Hy: | 
177 


relationship B=F(H) is determined by the relationship between Hg and Hy, Two cases are possible. 
1. Hg < Hg. This case corresponds to 


(20) 
BA, 


If Hs < H < Hg, there are no saturated portions of the core, and 


Ro 
| (kh —H,.) dR =p(H —H,). (21) 
I 
For Hy < H < Hy, where 


R alne 
Hy = = (He (22) 


is the field strength corresponding to complete core polarity reversal, it is easy to find the average core induction: 


B 


a—i . 
«ina 
Thus, if Condition (20) holds, the ascending arm of the hysteresis loop consists of a linear (central) portion 
(21) and two nonlinear portions, (17) and (23), at the beginning and end of the arm, 


2. Hs > Hg. ‘In this case 


(24) 


For Hy <H < Hs, the core has a portion with positive saturation (B,), adjacent to the inner perimeter, a 
portion with negative saturation (— Bs), adjacent to the outer perimeter and an intermediate, nonsaturated 
portion, Therefore, the average core induction is determined from the equation: 


where 

(#,+—*)ina (#,—-—*) ina 

By integrating, and then substituting the values for Ry and Rg, we obtain 
147 

B, = — (25) 


For Hy <H < Hy, the core induction Bg is determined from Formula (23), 


Thus, if Condition (24) holds, the central portion of the ascending (and the descending) arm of the hysteresis 
loop is characterized by a linear dependence of the induction on the field strength (25), However, unlike the 
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previous case (20), the differential core permeability on this portion, 


B 

(26) 
core” dil In o ‘ B, 
ull, 


does not equal the differential material permeability, and may be significantly less, 


It is interesting to note that, for pH¢ » B, and, in particular, for y—> oo, which corresponds to the ideal 
rectangular hysteresis loop, we have 


and Equations (25) and (26) take the following forms: 


B, = B,(— (25") 
B, 
core = . (26") 
With this, H, = Hy aad Hy = Hy, f.e., the ascending and descending arms of the core's hysteresis loop consist 
of just one linear portion (25°). 
If we set By =0 in Formulas (25) and (25"), and compute the corresponding values of H, we obtain the 
magnitude of the core’s coercive force Ho, which turns out to be larger than the material's coercive force. 


For example, from Formula (25'), we have 


a+i)ilna 
H.. 


a-4 


ti 


Fig. 8 


If the average core field strength is computed by Formula (5), then its coercive force turns out to be: 
Ho =Ho. 


14 2s 
GA, B, B 
In (1 + 2p) = 
uff, 
on; 
& gauss 
i 
1 
ly 
| 
/ 
-s000}' 
y 
16000 J | 
Fig. 7 
779 


Figures 7 and 8 give the static hysteresis loops for toroidal cores with a = 4, computed by the formulas 
given above, The dotted lines represent the material hysteresis loops, The hysteresis loop of Fig. 7 corresponds 
to the first.case considered above (B, = 15,000 gauss, H, = 0.2 oersted and 4 =10° gauss/oersted), and the loop 
of Fig. 8 corresponds to the second case (Bs = 15,000 gauss, He = 0.2 oersted and p = 10° gauss/ oersted), 


It is obvious from Fig. 7 that, if Condition (20) holds, the 
+4, Bt Hy) relationship B =F (H), on the ascending and descending arms of 
the hysteresis loop for a core with a sizeable value of a =Ro/Ry, 
a i ‘py differs significantly from a linear one, For example, for a =4, 
0 


the induction region,where the relationship B =F (H) is linear, 
consists of only 21% of the entire range of variation of the in- 
duction (2B,), 


The deviation of the relationship B = F (H) from the linear 
engenders the same type of infringement of “input-output” 
characteristic linearity as certain types of magnetic amplifiers 
with positive feed-back [1, 2]. 


“& It is clear from Fig. 8 that, if Condition (24) holds, an 
Fig. 9 increase in a shows its effects principally in a decrease in the 
differential permeability dB/dH on the “vertical” portions of 
the hysteresis loop of the core, which leads to a decrease in the gain and quality in the corresponding types of 
magnetic amplifiers with positive feed-back. 


By substituting in Formula (26) the values, typical for alloy 50NP, of B, = 15,000 gauss, H, = 0.2 oersted 
and » = 10° gauss /oersted, we obtain the following values of 1 core: 


a i—1.16 1.4 2 | 4 


10® | 4.5 10° |2.16 105) 1.1 10 


Yoore, 
gauss /oerst,} 
The formulas given above can also be used to compute the rectangularity coefficient of a hysteresis loop, 


which is one of the important parameters characterizing the quality of cores used in magneto-static triggers and 
memory devices, 


Re We assume that a field strength of Hp is required to re- 
7 verse a core's magnetic polarity from the state — B, to the state 
\ B,. Then, one ordinarily understands by the “rectangularity co- 
\ \ efficient" the ratio (Fig. 9): 
3 2 
as \ 
ar 
re The closer Re is to unity, the more reliably does the core 
a function in coincidence circuits and matrix memory devices, 
N\ p where a core must reverse polarity only if signals are applied 
simultaneously to two windings, and must not reverse polarity if 
a signal is applied to only one of the windings (see, for example, 
a ans : ay ame [2], paragraphs 13-18), A decrease of Re leads to an increase of 
the noise level and could lead to incorrect functioning of the 
Fig. 10 devices mentioned, 


As applied to our idealized hysteresis loop (Figs. 6-8), Hm = Hy, B(Hm) =B, and 


(27) 


f 


Figure 10 shows the dependence of Re on a, computed by Formulas (17) - (25) and (27), for three values 
of #H,/Bs. The curve (1 H_/Bs = 3) corresponds to an alloy of type 50NP and molybdenum Perminvar, with the 
values H, = 0,1 oersted, Bs = 15,000 gauss and p = 450,000 gauss/oersted, Even for =1.5, the rectangularity 
coefficient falls to 0.56. 


An increase of H, and a decrease of B, (for example, by the use of ferrite cores) allows Re to be increased 
for a given value of a . For #H_/B, =15 (curve 2), the rectangularity coefficient is close to the limiting value 
of Re, corresponding to # =0o, and any further increase 
8, of the ratio, » H¢/Bs, leads only to a small increase of 
Re, fora given a . For mp =a, we have Re =1 if 


a = 2, andRe= 


for a > 2 (curve 3), 

It should be expected, from a consideration of 
Figs. 1, 2 and 3, that the influence of the magnitude 
of a on the form of the limiting hysteresis loop of 
materials with a small ratio of residual induction B,, 
to saturation induction B,, will be significantly less 
than for materials with rectangular hysteresis loops, 


On Fig, 11, the solid curve is the descending 
arm of a limiting hysteresis loop, obtained experi- 
mentally for a core of alloy 79NMA, with a value of 
y a =1,.1, The dotted curve was obtained by a graphical 

oerst, 
| ap. integration of Equation (1) for the value a = 4, 


It is of interest to note that if, in Formula (26), 
- the value of B, is substituted for that of B,, this formula 


gives a value, Ut core = 146,000 gauss /oersted (for B, = 
= 5000 gauss, H, = 0,05 oersted and = 10° gauss /oersted), 
which practically coincides with the value of 1 core 
obtained above. 


5000 SUMMARY 
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Fig. 11 


The magnetic characteristics of toroidal cores, 
taking into account the influence of nonuniformity of 
magnetization along the cross section, can be obtained 
in the form of a mathematical relationship, B = F (H), 
if one has an analytic expression for the corresponding material magnetic characteristics, or may be obtained 
graphically, if the material characteristics are given in the form of experimental curves, 


The ratio of the outer and inner diameters of toroidal cores has an essential effect on their magnetic 
properties particularly when magnetic materials are employed which have rectangular hysteresis loops, or sharp 
corners in the basic magnetization curves, 


With an increase in the ratio of outer and inner diameters, there is observed a flattening of the corners 
on the magnetization curve, and a decrease of magnetic permeability in this region, With this, there is a de~ 
crease of the differential magnetic permeability on the “vertical” portions of the hysteresis loop, a departure 
from linearity of the relationship B=F(H) on these portions, and a decrease in the core’s rectangularity co- 
efficient. 


The variation of the magnetic characteristics of toroidal cores with increases of the ratio of outer and 
inner diameters can have an essential influence on the properties of magnetic amplifiers and other contact-less 
magnetic elements of automatic equipment, and must be taken into account in choosing the geometric dimen- 
sions of toroidal cores, 
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SOME PROBLEMS IN THE DESIGN OF AN ASYNCHRONOUS CLUTCH 
WITH A MONOLITHIC ROTOR 


T. A. Glazenko 
(Leningrad) 


Formulas for the maximum torque and critical slip are deduced, The effect of various 
design factors on the torque and the shape of the mechanical characteristic is determined. 
Relationships are given for the determination of the optimum number of pole- pairs for the 
asynchronous clutch, 


1. Introduction 


Asynchronous slip clutches (AC) with uniform rotors find application in the automatic control and re- 
gulation of the speed of electric drives for machine tools and transport machines, and also as actuating mechanisms 
in tracking systems, The use of an AC with a monolithic rotor permits the utilization of asynchronous motors 
with a short-circuited rotor in 3-4 kw controlled installations with heavy starting conditions, shock loading, and 
a repetitive short duration mode of operation [1]. 


The magnetic field in the AC is set up by the dc winding situated on the inductor (Fig. 1), The armature 
of the clutch consists of a cylinder made of a ferromagnetic material, 


The magnetizing force of the clutch Iw is constant and is independentof the relative rotor speed, Con- 
sequently, in setting up the electromagnetic field equations, the equivalent linear loading of the inductor A 
is taken as a basic quantity, and the field in the armature and gap of the AC is examined, 


Equations for the clutch's maximum torque M,,, which are obtained from the solution of the system of 
equations for the electromagnetic field in the gap and rotor of the clutch, are given in this paper, 


The derived relationships make it possible to evaluate the effect on the torque and the form of the 
mechanical characteristic of the AC, of various design factors: the dimensions of the magnetic circuit the 
outside diameter of the inductor D, the active length of pole | and air-gap § cm, the relative magnetic per- 
meability » h/cm, conductivity of the rotor steel y ohm=! cm!, and the number of pole-pairs in the inductor 


The inductor can be nonphaneropolar, phaneropolar, toothed with intersecting poles, and finger- toothed 
(Fig. 1). 


An attempt is made in the paper to compare the operating properties of the AC with various types of in- 

ductors, 

2. Clutch With a Nonphaneropolar Inductor 
4, Linear loading of the inductor, 


An inductor with not clearly defined poles makes it possible to obtain a close approximation to a sine 
wave in the curve of the distribution of the linear load A, and the normal component of induction By, in the 
gp, for ideal lost motion, 
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The ampere-turns of one pair of poles of such an inductor, expressed in amps, are found from the formula 


24 ma D 
0 


where x is the running coordinate on the evolute of the inductor's circumference, 


p 


Fig. 1. Mechanical schematics of asynchronous clutches, 

a) Clutch with a nonphaneropolar inductor; b) with a phanero- 
polar inductor: 1) inductor, 2) rotor, 3) field winding, 4) 
gap, 5) pole; c) with a toothed inductor and alternating 
poles; d) with a toothed inductor. 


Taking into account that the pole separation r = 1 D/2p, we obtain 
Iwp = AmaxD / p- 


Af we assume that the magnetic circuit of the clutch is unsaturated, then the amplitude of the nonlinear 
inductor load, expressed in amp/cm, can be expressed in the form 


2 
~ 
b 
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Here By max is the amplitude of the inductor in the gap, for ideal idle motion, when the angular fre- 
quency of the armature current we,=0, 6¢ is the size of the equivalent air gap in the clutch, determined 
during the design of the magnetic circuit, py =4 10°* h/cm is the magnetic permeability of air. 


In the design of a clutch, the amplitude of the induction in the gap, for ideal idle motion, is first picked, 
By max = 8000 ~ 12,000 gauss, and the necessary ampere-turns Iw, are determined, 


-b, Torque and optimum number of pole- pairs. 


The equations for the electromagnetic field in the gap and rotor of an AC with not clearly defined poles 
are set up under the usual assumptions [2, 3], 


™, kgm The joint solution of Maxwell's equations yields the 
following expression for the torque on the armature due to 
#4 the field: 
7 (2-+-V2)M 
ui—f V wet 


/ 
12 


where w is the relative angular speed of the rotor: 


o, .—@ 
as The relative angular speed for which the clutch 
torque reaches a maximum is: 
Pina” * rpm 
P 
400 800 1200 1600 @cr= Ong Se, = (3) 


Fig. 2, Calculated mechanical characteristics 
for a nonphaneropolar AC with D = 13,4 cm, 
=3,0 cm, 6 =0.03 cm, Iw, = 1500 amp, 


Here w,, and wa, are the angular speeds of rotation 
of the inductor and armature, s,, is the critical slip, 


=0.034cm, =10°, =0,312 x The maximum (critical) torque is 
wc 
x 10° ohm™! cm™4, Dl 
Mer = 0.266xp, A (4) 


where Ap is a coefficient which takes into account the variation of the magnetic permeability in the direction 
of propagation of the electromagnetic wave in the rotor [4]; Ap w 1.3 1.4, 


Formula (2) can be used for a qualitative evaluation and comparison of the shapes of the mechanical 
characteristics of different AC, 


A more accurate plot of the curve M = f (s) for an AC with a toothed inductor, within the range of 
operating speeds, can be obtained if one makes use of the formula given in [1]. 


Equation (2) indicates the presence of a maximum in the clutch’s mechanical characteristic, which 
corresponds to physical phenomena and is confirmed experimentally (Figs. 2 and 4), In fact, for synchronous 
rotation of the inductor and armature (the mode for ideal idle motion We) = Wjnq sp = 0) eddy currents are not 
induced in the armature and the torque is equal to zero, However, if w.;— then the depth of penetration 


1 
of the wave into the mass of the rotor z= —————————._ vanishes ther with the rotor losses and the torque, 
Gina toge 


Wind 
The value of the relative speed wjng8,, at which the maximum torque is reached depends on the dimen- 
sions of the clutch components (D, §) and the electromagnetic characteristics (u, y) of the armature material, 
The value of w,, increases sharply with increase in the number of pole~pairs or the length of the air gap, as 


well as with decrease of the clutch diameter, An increase in w,, with M., = const corresponds to a smoothing 
of the AC mechanical characteristic, 


It should be noted that in clutches with a sharply nonsinusoidal distribution of the equivalent, linear load 
on the inductor in the gap, or with an appreciable number of pole-pairs, the maximum torque is usually reached 
for amounts of slip lying well outside the limits of the operating range of changes in the armature speed. 


The expression for the critical torque follows from Relationships (1) and (4): 
Mex = 10.5 x 10-* p(Iw,)* Ay. (5) 


If the dimensions, ampere-turns Iw,, and relative speed wjnq8, for which the nominal torque is reached, 
are choten, then the most profitable number of pole -pairs Pope can be determined, 


An increase in the number of pole~pairs is accompanied by an increase in the electric speed of eddy 
currents in the rotor w.) = WingSp. At the same time the armature heat losses increase, while the clutch torque 
increases up to a definite limit. With further increase in p the torque begins to decrease, i.e, M—> 0 as 
p> o, due to the ever growing effect of the armature reaction and the decrease in the depth of penetration 
of the electromagnetic wave into the rotor. 


To determine the optimum number of pole- pairs it is convenient to transform Expressions (2), (3), and 
(5) to the form 


M =3.58 x 10 


= Cf (K » (6) 


where 


The torque resches « maximum when =0, o 2.5 Kp V2p=0. 


In this way, the optimum number of pole-pairs depends on the value of K which, in its turn, is determined 
by the clutch parameters, 


The relative speed wings in (7) is chosen from the conditions for a stable transmission of torque, Usually 
one takes Wing (1.5 3) wing Sp- 


The curve Pope = f (k) (Fig. 6, curve a/D-> 0) is used to determine the optimum number of pole- pairs 
for an inductor with not clearly defined poles, 


38. Clutch With Phaneropolar and Toothed Inductor 


a, Linear load on inductor, 


For the case when we} =0, the distribution of induction in the gap, for a phaneropolar, toothed inductor 
with intersecting poles, has a trapezoidal form, The corresponding curve, of the distribution of the equivalent 
linear load A= f (x) along the circumference of the inductor is given in Fig. 3a. 


The maximum value of the linear load on the inductor is A‘,,,, = 2B) max 5e/Ho@, Where a is the cal- 
culated distance between the poles of the inductor. For a clutch with a phaneropolar inductor 


Iw 
Ainax = p/w (8) 


P 


~ 
K , p? K , p® 
p* 
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For a clutch with a toothed inductor (with intersecting poles) 


Iw 
A'max= = Iw, (9) 


where Iw; is the total ampere-turns packed on the inductor, a is the distance between the intersecting teeth. 
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Fig. 3. Distribution of the equivalent linear load on the inductor and the normal 
component of induction when w=0: a) for a phaneropolar and toothed (with alter- 


nating poles) inductor, b) for a toothed (finger-type) inductor, 1) Armature; 
2) coil; 3) inductor, 


To obtain a certain value of induction for idle motion, the clutch with toothed inductor must contain p 
times less ampere-turns than a clutch with clearly defined poles, 


For a toothed (finger-type) inductor, the induction changes from By max t© Bymin over the length of a 
very small portion of the gap's circumference (Fig. 3b), Here a =0,1 - 0.3 cm, usually, 


In this case the maximum value of the equivalent linear load can be found from the equation 


x a 


ma 


Iw tot 
where Iw = 


The active length of axis in the finger-type of inductor is? = 2, (ig. 1). 
ng P 


la, B. Rozman [1] proposed a diagonal displacement of poles in the finger-type of inductor, Such an 
arrangement weakens the edge effect on the inner sides of the poles, which results in an increased torque, 


The distribution of the equivalent, linear load along the circumference of a toothed and phaneropolar 
inductor can be represented by a trigonometric series of the form 
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where k =1, 3, 5,... — is the order number of the harmonic of A, Ay max is the amplitude of the k-th harmonic 


component of the curve A= f (x): 


ma a 
nD 


T and number of 


for the clutch, 


For a clutch with a phaneropolar and toothed inductor, the torque is determined by the action of all the 


harmonic components of the linear load, 


M kgm 
20 
16 
p=2 


j@0 [200 th 


Fig, 4, Calculated and experimental character- 
istics of an AC with a phaneropolar inductor 
with D=13,4 cm, 6 =0,04sec, 1 =3,.0 cm, 
= 0.31210" ohm™! =10°, 5, 
= 0.071 cm, Iwp =1500 amp. The circles 
denote experimental points, 


a 


a5 shop 


Fig. 5, Calculated mechanical characteristics 
of an AC with a toothed (finger-type) inductor, 
for D=14,8 cm, I = 21 6 =0,03 cm, 
Yeo =0.312- 10° =10°, 

5 ¢ = 0.034 cm, P= Popt = 6, = 800 turns. 


The torque's component determined by the action 
of the k-th harmonic of A is equal to 


cr k 
The angular velocity, for which the component 


of torque due to the action of the k-th harmonic of A 
reaches a maximum, is equal to 


= (12) 


16 (pk)*3? 


The maximum torque component due to the 


_ action of the k-th harmonic component of A is equal 


to 
Meg x = 10.5 x 10°F (14) 


Exp..:ssions (13) and (14) show that M,,;,, due 
to the action of the highest harmonics of A, even for 
a linear load to Amaya,» is k times 
smaller than the torque due to the action of the first 
harmonic component, Slip for which M,,, is reached, 
is k* times greater than Usually, scrg, Scrg» 
lie far beyond the limits of the values of slip which 
occur in the clutch under operating conditions, Because 
of this, in the range of operating speeds, the torque is 
mainly due to the action of the first harmonic com- 
ponent of A, 


The mechanical characteristics can be obtained 
by summing the ordinates of the mechanical character- 
istics My, = f (w) which are plotted for each harmonic 


of A (Figs. 4 and 5), M,= >) Ms. Some diver- 
k=1,3... 

gence between the calculated and experimentally de- 

termined characteristics (Fig. 4) is explained by the 

fact that in the derivation, variation of the relative 

magnetic permeability of the rotor with we] was not 

taken into account, 


ue It is a characteristic of this type of clutch that the maximum in its mechanical characteristic is attained 


at considerably higher relative velocities, and is strongly smoothed out due to the action of the highest harmonics 


of A (Figs.4 and 5), In accordance with (4) and (11), the critical torque, which is determined by the fundamental 
harmonic of A, is equal to 


M,,.1=17 x = ci(p, 
| 


In this way, the smaller the ratio — is, the larger is Mor, and, in the limit, 


47 10° p (Imp)? An white 
To increase the torque, the spacing between 
Pa poles (or intersecting teeth) must be minimized, 
a0 The minimum value of a is limited by the increase 
4% as of the leakage flux between the poles, Usually one 
2 ee selects a =(10-20)8. This condition is always met | 
ge ra in clutches with finger-type inductors, since in these | 
p= ~ ath the change in induction in the gap occurs over a . 
es =aze5' | very small portion of the length of the inductor, The 
optimum number of pole-pairs p,,, can be deter- 
fone: the condision fer tesque by 
equating to zero the derivative, with respect to p, 
Fig. 6, Curves of the relationship Pop; = fh 
with a/D =const, 
sin? 
M =56.3x10° 2 ( =c''(4,X). 


The optimum number of pole- pairs depends, in this case, not only on the magnitude of K, determined by 
Equation (7), but also on the ratio 4 (Fig. 6). 


e The curves Pop: = f (k), for several values of > + were plotted to determine popt (Fig. 6). 


Example, Let us determine the pop; for an AC of the PMS-0.4 type, with a toothed inductor, 
Given: D=14,4cm, =0,03 cm, a=0.3.cm, nind = 1340 rpm, =0.1, =0.312-10° cm", 
10°, h/cm. The calculations are carried out in the following order: 


= 2 ing? n= 256 /secm!  K = 48: = 5.4, 
/ D =0,3/ 14,4 = 0.0208. 


From the curve Pop, = f () (Fig. 6) corresponding to = = 0.0208, for the value K = 5,4, we find 
Popt = 8-7. 
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| 
| 
| 
7189 


SUMMARY 


The maximum torque in an AC with a phaneropolar inductor with intersecting teeth is greater than in 
clutches with finger-type and nonphaneropolar inductors. 


It is necessary to select the smallest ratio of ~ to increase the torque in an AC, 


The mechanical characteristics of nonphaneropolar clutches have very sharply defined maxima. 


The sharp mechanical characteristics of an AC with a phaneropolar and toothed (with intersecting poles) 
inductor increase its efficiency of transmitting the nominal torque. A clutch of this type is most useful for 
applications where large values of transmitted torque are required, 

In clutches with diameters D <(10 - 20) cm, the distribution of the required ampere-turns of excitation, 
for a phaneropolar and toothed inductor with intersecting poles, leads to complications and often makes it 
necessary to increase the clutch dimensions, In this case, a clutch of very simple design, with a finger-type 
inductor, possessing relatively smooth mechanical characteristics (Fig. 5), should be used, 


It is also expedient to use an AC with a finger-type inductor in installations with a wide range of speed 
variations, where the clutch dimensions are determined not by the value of the transmitted torque, but by the 
necessary area of the heat-emitting surface. 


When used in systems for the automatic control of the speed of rotation of an actuating device by changes 
in the driving current, the AC must contain a speed feedback loop, 


LITERATURE CITED 
(1) Ia, B, Rozman, “Electromagnetic slip clutch for the stepless speed control of an actuator,” [In 
Russian) (Institute of Techno-economic Information, AN SSSR, 1954), 


_ [2] E, M, Lopukhina, “Study of an asynchronous motor with a hollow, cylindrical rotor," Electrichestvo 
(USSR) 5 (1950). 


[3] Alfred S,Cutman, “Designing an eddy current clutch servo,” Mach, Design 10 (1950), 


[4] L, R, Neiman, in the book: Surface Effect in Ferromagnetic Objects [In Russian] (Gosenergoizdat, 
1949). 


Received July 27, 1957 


COORDINATION CONFERENCE ON AN AUTOMATED ELECTRICAL 
AC ACTUATOR WITH A SELF-CONTAINED SUPPLY 


A coordination conference, organized by the Power Institute, A. N, Azerb, SSR and the Institute of Auto- 
mation and Remote Control, AN SSSR, was held on October 7-9, 1957, in Baku, concerning an automatic, 
electrical, alternating current actuator with a self-contained supply. 


The conference was attended by 120 delegates from 44 organizations in 14 towns in the Soviet Union, 
including representatives of the IAT AN SSSR, the Electromechanical Institute of AN SSSR, the Moscow Power 
Institute, the Power Institutes of AN Uzbek SSR, AN Laty. SSR, the L'vov Polytechnical, Azerbaidzhan Poly- 
technical, Azerbaidzhan Industrial, Cheliabin Polytechnical, Khar'kov Polytechnical, and Erevan Polytechnical 
Institutes, the Orden Lenina Military Academy of Tanks Corps, Giproneftmasha, Ts, K, B, “Electroprivod,” 
Bashenergonefti, the Ministries of Petroleum Industries of Azerb, SSR, etc, 


Thirty papers were delivered at the conference, The proceedings of the conference were divided into two 
parts: 


a) self-contained, alternating current actuators; 
b) automatic regulation and control of alternating current actuators, 


In the introductory speech, the vice-president of A, N, Azerb SSR, Professor Z, 1. Khalilov noted the 
significance of the conference, and formulated the problems which have to be solved. 


In the paper of Academician V, S, Kulebakin, "Problems of automatic, combined regulation of electrical 
actuators,” the theory and methods of designing systems with combined regulation by means of reaction to de- 
viations and changes, were elucidated, It was shown that by using the combined method of control and regulae- 


tion, separate coordinates of an automatically- operating system can be made almost independent of external 
actions or loads. 


The theory and practice of throttle and frequency control systems was the subject of the second paper of 
Academician V, S, Kulebakin, “On the application of saturable reactors and rectifiers in widely controlled 
actuator systems." The speaker discussed the main circuits of direct current actuators, supplied from single- 
phase and three-phase power supplies, and controlled by means of saturable reactors, static characteristics of 
electric actuators, as well as circuits with saturable reactors and rectifiers for the frequency control of alternating 
current motors, 


The results of scientific studies of self-contained, alternating current actuators carried on at the Power 
Institute of A, N, Azerb, SSR were presented by the Corresponding Member of A, N. Azerb. SSR, A. A. Efendizade, 


The Institute has carried out studies of the automation of drill- feed on petroleum installations using rotary 
boring, and the automation of the brake machine control in drilling oil bores, A study was made of the opera- 
tion of an asynchronous motor when supplied by a synchronous generator of commensurate power, as well as & 
system for the regulation of the excitation of synchronous generators in self-contained installations, 


The paper of T. Z. Portnoi (Ts, K, B, "Elektroprivod") was devoted to self-contained electrical drives for 
machines in drilling installations, He presented the results of industrial tests of a direct current, diezel-electrical 
drive unit, compared direct current and alternating current diesel electric actuators, and formulated the problems 
in designing diesel-electric drives for drilling installations, 
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The report of V, N, Bogoiavlenskii (IAT AN SSSR) “Regulated asynchronous actuator, using planetary 
transmissions, in self-contained supply systems" was devoted to a frequency control system, incorporating 
electromechanical transmission of alternating current, developed for an electric tractor, 


The report of V. G, Vovka (Construction Bureau VNII) was also devoted to the application of electro- 
mechanical transmission, 


Supply systems with frequency control by means of commutating generators was the subject of the papers 
read by representatives of IEM AN SSSR, 


The report of V, P, Andreev and V, A, Prozorov, “An automated, alternating current actuator with very 
smooth starting and stopping,” treated the problems of automating the control of an electric actuator with large 
flywheel loads, by means of an alternating current motor- generator system using a compensated commutator 
generator, 


The reports of V. V, Rudakoy “The work of EM AN SSSR on quick-response, electromechanical fre- 
quency converters,” and E, M, Smirnov “A comparison analysis of some electromechanical systems for the 
excitation of a compensated, commutator generator,” examined various electromechanical frequency converters 
with an inverted or normal commutator, supplied from the windings of electrical machines or static transformers, 


In the report of A, A, Dartau, "On the selection of an optimum system for the excitation of a commutator 
generator,” experiments carried out on resonant systems for the excitation of commutator generators were described, 
An optimum type of commutator generator with resonance switching of field coils, was recommended, 


lu. M, Aleksandrov, in his paper “Single-phase, regulated, alternating current actuator with frequency 
control,” described investigations of an ac motor-generator set suitable for mechanisms not subjected to sharp 
load surges, and consisting of a single-phase commutator generator, based on an amplidyne, and an asynchronous 
motor connected single-phase, 


The report of V, M, Mamedov, “The operation of a longitudinal field, dynamoelectric amplifier with ac 
self-excitation,® set forth the results of work directed at the design of alternating current amplifiers based on 
standard types of direct current, dynamoelectric amplifiers. 


S. V. Strakhov (MEI) delivered a report on the calculation of the static stability of diesel- electrical, 
alternating current actuators, A method is given for the determination of the static stability in the absence, as 
well as in the presence of voltage and current feedback loops. It is shown that the limit of the static stability 
for an asynchronous motor supplied by a generator of commensurate power- rating, shifts appreciably in the 
direction of increasing slip, in comparison to the case of a motor fed by an infinite power network, 


E, A, lakubaitis (Power and Electrica] Engineering Institute A. N, Latv. SSR), in his report, "Self-contained 
power station for a railway passenger coach," presented results of the operation of systems for the supply of 
electricity to passenger coaches, in the form of rectified ac, using a synchronous, contactless generator, and a 
system for the automatic control of the excitation of a self-exciting synchronous generator with variable rpm, 


The problems of stabilizing the voltage of a synchronous generator, operating in a self-contained system, 
were treated in the reports of M, M, Rasulov (ENIN A, N, Azerb, SSR) “A system for the regulation of the ex- 
citation of synchronous generators in a diesel-electric actuator,” V. A, Glebov (RIL Zh, T.) “Forced starting 
of asynchronous motors in self-contained electrical systems," V. T. Zagorskii (L'vov Polytechnical Institute) 
**Capacitor start of short-circuited asynchronous motors from generators of commensurate power." 


The report of M, M, Rasulov formulated the requirements, based on theoretical and experimental investi- 
gations, imposed on regulators of the excitation of synchronous generators in self-contained installations, The 
author presented the results of investigations of new circuits using saturable reactors for the control of the ex- 
citation of a generator. 


V. A, Glebov's report presented the results of an investigation of the starting of a short-circuited asyn- 
. chronous motor from a synchronous generator of commensurate power, for pulse starting, starting through an 
SO auto-transformer, and through a resistance, 


V. T. Zagorskii described, in his report, the results of an analytical and experimental investigation of a 
synchronous generator — asynchronous motor system with series-connected, fixed capacitors, The speaker gave 
the conditions for the self-excitation and self-drive of the motor. 


In the report of M, A, Nabiev (ENIN AN Azerb SSSR), “Investigation of transient processes in a syn- 
chronous generator — asynchronous motor system with regulated excitations," the theory of the control of a syn- 
chronous generator's excitation system was presented. 


V. lu, Nevraev (IAT AN SSSR) in his report “Analysis of transient processes in, and choice of structure 
for a system for the automatic voltage control of a self-contained synchronous generator — asynchronous motor 
installation,” described a semi-graphical method of calculating the transient processes in an exciter— synchro- 
nous generator — asynchronous motor system, 


V. A, Parail (Odessa Polytechnical Institute), in his report “An alternating current generator- motor system 
for diesel-electric actuators,” presented a method for obtaining sharp regulation characteristics of an asynchro- 
nous motor, operating from a synchronous generator, by acting on the excitation circuit of the generator's exciter, 


New types of supplies for boring installations were described by the representatives of Giproneftemasha, 
L, I, Shuturman in “Drive of drilling machines from mobile. gas- turbogenerators,” and I, N, Sulkhanishvili in 
**Diesel-electric, alternating current, drill drive, type DEB- 2," 


In these reports the authors outlined the aspects of designing a gas, turboelectric, drilling unit, pointed 
out the basic requirements imposed on diesel-electrical, drilling actuators, and presented the basic facts about 
an experimental drilling installation with gas-turboelectrical units, and about the circuit of a DEB-2 diesel- 
electrical actuator, 


The report of I, I, Giultmamedov (ENIN AN Azerb SSR) “Mechanical transient processes in an asyn- 
chronous motor supplied by a self-contained synchronous generator," contained the results of investigations of 
mechanical, transient modes in an asynchronous motor, in a self-contained system, under various load conditions, . 
and the equations describing these modes were given. 


B, A, Listengarten, in his report “The geometrical location of currents and voltages in an asynchronous 
motor operating from a synchronous generator of commensurate power," determined the location of currents and 
voltages in an asynchronous motor, taking into account the motor’s magnetizing circuit, when working from a 
phaneropolar, constant excitation, synchronous generator, 


In the second section, a report was given by Ia, B, Kadymov (ENIN AN Azerb, SSSR) “Stability criteria 4 
for automatic control systems, with distributed parameters,” devoted to the popularization of frequency methods 
in the analysis of the stability of systems with distributed parameters, The report outlined an engineering . 
method of investigating the stability of systems, described by complicated transcendento- algebraic equations . 
(for example, systems with two delays). | 


Several reports were devoted to electrical actuators applicable to the petroleum industry, 


K, N, Kulizade (Azerb. Industrial Institute), in his report "The application of low power, synchronous 
motors in installations of the oil industry," presented the results of the experimental application of a synchronous 
motor, with mixed excitation, and metallic rectifiers as a drive for pumping rigs on one of the oil fields in 
Azerbaidzhan, 


The report of In, S, Kengerlinskii (ENIN AN Azerb, SSR) was devoted to the investigation of an automatic 
control system for tool-feed, as well as the results of investigations of a power regulating circuit, and a set of 
requirements, imposed on systems used in rotary drilling, were established. 


A. M, Omarova (ENIN AN Azerb, SSR) described in her report the results of investigation and industrial 
trial of a brake machine (a synchronous generator operating in the region of dynamic braking) with a system for 
the automatic control of the lowering of an instrument (continuous and pulse), 


The report of A, A. Mezhlumova (Azerb, Polytechnical Institute) “Electric drill — an electric drilling 
machine perfecting the technique of drilling,” was devoted to the description of a machine which combined in 
one the electrical supply for the rotary table and the rotary table itself. 


D. P, Petelin (IAT AN SSSR) delivered the report “Investigation of a system for the automatic control 
of the frequency and voltage of a synchronous generator" in which an investigation is made, based on the theory 
of multiple-connection control, of a system for the automatic control of the frequency and voltage-of a syn- 
chronous generator, under the condition of any stable accuracy of each of the controlled quantities, 
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Arrangements for the oscillographic display of the angle @ in transient modes of synchronous machinery 
were described in the reports of G, M, Safarov (Azerb. Industrial Institute) and M, la, Ginzburg (ENIN AN 
Azerb. SSR). 


Thirty-two persons took part in the discussion of the reports and the conclusions of the conference. 


Several measures were planned in the conclusions of the conference, of which the following were the 
most important, 


1, To intensify theoretical work on an automated electrical actuator with a self-contained supply, in 
scientific institutions of the Academy of Sciences, USSR, the Academies of Sciences of Soviet Republics, as 
well as in schools of higher education and branch scientific institutes, 


2. To consider one of the basic problems in the field of self-contained diesel-electric drives to be the 
investigation of special versions of motors and generators which would ascertain the proper operation of actuators, 
and the investigation of complex systems for the control of the thermal and electrical parts of these actuators, 
The case of the simultaneous operation of several self-contained generating units with a common, sharply 
varying load, should also be considered, 


3. To consider one of the most important theoretical and practical problems in the field of self-contained 
diesel-electrical actuators, to be the design of automatic systems for the remote control of prime movers; as 
well as the determination of optimum conditions and requirements for such control, 


4, To devote proper attention, in the design of diesel-electric actuators, to the development of systems 
for the control of electric motors, making wide use of contactless automation, permitting the shaping of optimum 
load diagrams for the generator and prime mover. 


5. In the investigation of transient processes in self-contained systems, it is necessary to take into account 
changes in the angular speed of the synchronous generator's rotor, and not to limit oneself to the development of 
semigraphical methods of investigation. Computers should be used widely. 


la, B, Kadymov 


ERRATA 


In the note of M, A, Aizerman and F, R, Gantmakher, “On the stability of periodic modes in nonlinear 
systems with piecewise-linear characteristics," Automation and Remote Control 19,6 (1958), 


1, On p. 594 replace (1) everywhere by [4]. 


2. Number [4] in the Literature Cited should read: 
ed 


[4] M, A, Aizerman and F, R, Gantmakher, “The stability of the linear approximation to the periodic 
solution of systems of differential equations with discontinuous right sides," Applied Mathematics and Mechanics 
(USSR) 21, 5 (1957). 
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